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Photodesorption  of  a molecule  from  a metal  surface,  following  absorption 
of  a visible  or  UV  light,  involves  electronic  transitions  of  the  desorbing  species 
coupled  to  the  lattice  vibrations  and  electron-hole  excitations  of  the  substrate.  This 
dissertation  presents  a general  treatment  of  those  phenomena,  based  on  a stochastic 
theory  of  localized  perturbations  in  an  extended  system.  The  treatment  is  later 
specialized  for  diatomic  adsorbates,  in  particular  carbon  monoxide  adsorbed  on 
the  (001)  surface  of  nickel. 

The  extended  molecular  system  Hamiltonian  is  divided  into  a term  for  the 
localized  region  of  primary  degrees  of  freedom  (DFs)  affected  by  the  desorption, 
coupled  to  the  secondary  DFs  of  the  substrate  that  acts  as  a time  evolving  bath. 
A self-consistent  treatment  gives  an  effective  Hamiltonian  for  the  primary  DFs 
that  accounts  for  energy  fluctuation  and  dissipation  in  terms  of  the  properties  of 


adsorbate  and  substrate.  The  assumption  of  Brownian  motion  in  the  substrate 
shows  the  physical  origin  of  dissipation  in  the  primary  effective  Hamiltonian  in 
terms  of  time  correlation  functions  (TCFs)  of  lattice  displacements  or  electronic 
density  fluctuations.  A diabatic  electronic  representation  is  used  to  eliminate 
momentum  couplings  between  adsorbate  electronic  states.  The  bath  dynamics  is 
studied  for  lattice  vibrations  and  electronic  excitations.  To  include  interaction 
with  light,  a radiation  field  term  and  a system-field  coupling  term  are  added  to 
the  Hamiltonian. 

A time  evolution  algorithm  is  developed  for  a system  with  two  DFs  and 
two  coupled  complex  potential  energy  surfaces  (PESs).  and  applied  to  the 
CO(ads)Ni(OOi).  The  time  dependent  Schrddinger  equation  for  the  primary  DFs  is 
solved  numerically,  using  time  evolution  of  wavepackets.  The  primary  DF  wave- 
function  and  the  transition  dipole  amplitude  arc  propagated  over  time  to  obtain  the 
autocorrelation  function  and  the  dipole-dipole  TCF,  respectively;  physical  initial 
conditions  are  specified  and  discussed  for  both  cases.  The  lifetimes  of  desorbing 
species'  quantum  states  are  obtained  from  the  autocorrelation  functions,  while  a 
Fourier  transform  of  the  dipole  TCF  yields  a generalized  desorption  lineshape  for 
non-Hermitian  dissipative  Hamiltonians. 

The  primary  PESs  and  their  couplings  were  parametrized  semicmpirically, 
using  experimental  data  and  the  results  of  the  INDO  electronic  structure  calcu- 
lations. Widths  of  the  PESs  and  couplings,  represented  by  imaginaty  potential 
functions,  were  estimated  from  the  data  on  desorption  lifetimes. 


The  results  of  the  lime  propagation  show  that  the  energy  dissipation  in  the 
substrate  plays  a crucial  role  in  the  dynamics  of  the  desorbing  species:  both  at 
short  limes  when  it  influences  the  transitions  between  the  primary  electronic  states, 
and  at  late  times  when  it  determines  the  photodesorption  lifetime.  The  desorption 
lineshape  results  show  a peak  positioned  in  qualitative  agreement  with  a previous 
experimental  result;  there  is  also  additional  structure  attributable  to  the  internal 
vibrational  excitations  of  the  desorbing  species. 


CHAPTER  1 
INTRODUCTION 


Statement  of  Ihe  Problem 

Photodesorption  of  carbon  monoxide  from  the  nickel  (001)  surface  is  studied 
in  this  dissertation  as  an  example  of  a broader  class  of  problems  including 
extended  molecular  systems  interacting  with  light.  These  extended  systems  may 
be  polyatomic  molecules  or  molecules  in  solutions  that  interact  strongly  with 
the  solvent;  this  work  focuses  on  small  (diatomic)  molecules  adsorbed  on  metal 
surfaces.  The  key  idea  of  its  approach  is  that  the  physical  processes  in  the 
adsorbate-substrate  complexes  arc  localized,  i.e.  describable  by  a small  number 
of  coordinates  or  degrees  of  freedom  (DFs)  of  the  adsorbate  molecule  interacting 
with  the  metal  substrate  that  contains  a very  large  number  of  DFs  describing 
different  excitations  in  it.  A physical  model  is  developed,  in  which  the  desorption 
follows  absorption  of  visible  or  UV  light  and  is  preceded  by  an  electronic  transition 
in  Ihe  adsorbate-substrate  complex,  coupled  to  the  lattice  vibrations  and  electron- 
hole  excitations  of  the  substrate.  One  can  thus  introduce  the  concept  of  desorption 
induced  by  an  electronic  transition.  The  localized  character  of  the  molecule-metal 
interface  processes  makes  it  possible  to  classify  the  DFs  in  the  system;  those  that 
describe  the  localized  processes  that  lead  to  the  desorption  constitute  the  primary 
DFs  of  the  system,  and  they  are  coupled  to  a very  large  number  of  secondary 
DFs  of  the  substrate. 


The  aim  of  this  study  is  to  solve  the  time  dependent  Schrodingcr  equation 
(TDSE)  for  the  adsorbate-substrate  complex-  The  concept  of  primary  and  sec- 
ondary DFs  is  used  to  partition  the  total  Hamiltonian  of  the  complex  accordingly. 
It  is  then  assumed  that  the  primary  and  secondary  DFs  are  coupled  in  a self- 
consistent  field  (SCF)  manner,  and  the  total  wavefunction  of  the  system  is  written 
as  a product  of  primary  and  secondary  wavefunctions  that  depend  upon  the  re- 
spective DFs  only.  This  formulation  gives  separate  TDSEs  for  the  two  sets  of  DFs 
and  this  work  focuses  on  solving  the  primary  TDSE  with  an  effective  Hamiltonian 
that  incorporates  the  effects  of  the  substrate.  The  substrate  dynamics  is  modeled 
stochastically  for'  different  types  of  excitations,  without  the  explicit  solution  of 
the  secondary  TDSE;  this  treatment  leads  to  energy  dissipation  in  the  primary 
DF  dynamics  that  is  included  as  a non-Hermitian  part  of  the  primary  effective 
Hamiltonian.  The  primary  TDSE  is  then  solved  by  means  of  wavcpackct  propa- 
gation; the  propagation  algorithm  was  constructed  to  include  the  presence  of  two 
electronic  states  of  the  adsorbate-substrate  complex,  and  the  effects  of  dissipation. 

Previous  Work  on  Dialomics  Adsorbed  on  Metal  Surfaces 

The  theoretical  and  experimental  study  of  desorption  induced  by  electronic 
transitions  (DIET)  has  been  a topic  of  considerable  interest  in  recent  years,  with 
several  literature  reviews  available  (Talk  et  al.  cds.  I9S3.  Brenig  and  Menzel 
eds.  1986,  Stulcn  and  Knotck  cds.  1988,  Betz  and  Varga,  cds.  1990,  Burns 
et  al.  cds.  1993].  In  particular,  the  stimulated  desorption  has  been  discussed 
in  reviews  of  experimental  work  (Lin  and  Ertl  1987,  King  and  Cavanagh  1989. 
Avouris  and  Walkup  1989],  Experimental  results  have  been  obtained  in  particular 


for  CO  adsorbed  on  Ni  using  the  electron  energy  loss  spectroscopy  (EELS)  and 
the  ultraviolet  photoelectron  spectroscopy  (UPS)  (Blyholder  1974,  Akimoto  ct  al. 
1979,  L6pez-Sancho  and  Rubio  1981). 

Recent  experimental  results  using  time-resolved  spectroscopy  have  probed 
the  substrate  dissipation  phenomena  in  mote  detail,  and  provide  new  information 
about  the  role  of  metal  phonon  and  electron-hole  excitations  in  systems  such  as 
CO  adsorbed  on  Cu  [Prybyla  et  al.  1990,  1992,  Cavanagh  et  al.  1993)  and  CO 
adsorbed  on  Pt  (Germcr  et  al.  1993  and  1994].  The  most  recent  theoretical  work 
[Newns  et  al.  1991,  Head-Gordon  and  Tully,  I and  n 1992]  deals  with  these 
time-resolved  experiments. 

The  basic  mechanism  of  stimulated  desorption  has  been  known  for  some  time, 
under  the  name  of  the  Mcnzcl-Gomcr-Redhead  (MGR)  model  [Menzcl  and  Gomer 
1964,  Redhead  1964).  This  model,  originally  constructed  for  atoms  adsorbed  on 
surfaces,  consists  of  two  electronic  slates  of  the  adsorbate-substrate  complex:  a 
ground  state  with  a potential  well  and  a repulsive  excited  state.  Absorption  of 
visible  or  UV  light,  or  scattering  of  e.g.  an  electron  beam,  causes  the  electronic 
transition  into  the  repulsive  excited  state  and  the  system  continues  to  evolve  over 

A similar,  but  more  complicated  desorption  mechanism  is  the  desorption 
induced  by  multiple  electronic  transitions  (DIMET).  It  occurs  when  the  excited 
electronic  slate  is  also  a bound  state,  so  the  diatomic  cannot  directly  desorb  in  it. 
Instead,  it  deexcitcs  back  into  the  ground  electronic  slate,  but  into  a distribution  of 
its  internal  (mostly  vibrational)  states.  This  process  can  be  repeated  several  times 


if  the  electronic  elicitation  is  continuously  induced  e.g.  by  light  pulses  of  visible 
or  UV  frequencies  corresponding  to  the  electronic  excitation  energies,  until  the 
diatomic  eventually  desorbs  in  its  ground  electronic  state  with  some  vibrational 
distribution.  Examples  of  this  behaviour  are  NO  adsorbed  on  Pt  [Gadauk  et  al. 
1990|  and  on  Pd  [Misewich  et  al.  1992). 

Notably,  the  MGR  model  docs  not  incorporate  the  continuum  properties  of 
the  substrate,  nor  does  it  include  in  its  original  version  the  internnl  (vibrational 
and  rotational)  primary  degrees  of  freedom.  In  more  recent  theoretical  studies, 
the  continuum  of  substrate  excitations  has  been  included  as  a finite  width  of  the 
primary  potential  energy  surfaces  (PESs).  The  width  is  expressed  as  an  imaginary 
potential  in  scattering  theory  [Micha  1976).  Its  form  has  been  derived  e.g.  from 
the  propagator  theory  [Brcnig  1976),  or  more  recently  in  theory  of  adsorbate 
spectroscopy  [Cederbaum  and  Tarantelli  1993),  It  can  be  shown  that  the  form 
of  these  imaginary  potentials  is  directly  related  to  the  stochastic  properties  of  the 
substrate  dynamics  [BefcSid  and  Micha,  to  be  publ.].  The  vibrational  motion  of 
diatomics  adsorbed  on  surfaces  has  also  been  amply  studied,  in  particular  for  CO 
adsorbed  on  Ni  [Richardson  and  Bradshaw  1979,  Black  et  al.  1981  and  1982, 
Parra  and  Micha  1986,  Micha  and  Parra  1988),  on  Cu  [Persson  and  Persson  I and 
II  1980,  Head-Gordon  and  Tully  1992)  and  on  Pt  [Persson  and  Rybcrg  1989). 


The  rest  of  this  dissertation  focuses  on  the  specific  problem  of  carbon  monox- 
ide adsorbed  on  the  Ni(00!)  surface. 


Chapter  2 contains  the  physical  model  for  this  extended  system,  based  on 
earlier  theoretical  [Sung  and  Hoffmann  1985,  Hoffmann  1988,  Drakovaand  Doyen 
1990)  and  experimental  [Akimoto  et  al.  1979,  Ldpez-Sancho  and  Rubio  1981) 
studies.  The  localized  character  of  the  photodesorption  motivates  the  concept  of 
partition  into  the  primary  and  the  secondary  region  [Kumamoto  and  Silbey  1982, 
Miller  1987,  Gerber  and  Rainer  1988,  Micha  and  Srivastava  1989,  Srivastava  and 
Micha  I and  II  1991,  Coalson  1989,  Messina  and  Coalson  1991],  and  the  effective 
Hamiltonians  in  terms  of  the  degrees  of  freedom  of  both  regions  are  derived  within 
the  self-consistent  field  treatment.  This  treatment  shows  the  connection  between 
the  energy  dissipation  in  the  secondary  region  as  a continuum  of  surface  states  and 
the  imaginary  potentials  introduced  phenomenologically  through  scattering  theory 
[Micha,  76).  Time  correlation  functions  are  introduced  at  this  point,  as  well  as 
their  SCF  formalism.  Physical  properties  such  as  lifetime  of  desorbing  states  and 
photodesorption  lineshape  arc  expressed  in  terms  of  the  TCFs.  and  the  lineshape 
is  derived  for  systems  with  energy  dissipation.  The  detailed  derivation  for  a 
transition  probability  with  a dissipative  Hamiltonian,  related  to  the  lineshape,  is 
presented  in  Appendix  A. 


The  stochastic  dynamics  of  the  secondary  DFs  is  studied  in  detail  in  Chap- 
ter 3,  for  the  cases  of  surface  phonons  and  electron-hole  excitations.  The  lattice 
vibrations  are  described  through  the  time  dependence  of  substrate  normal  mode 
coordinates  and  conjugate  momenta,  while  the  electron-hole  excitations  are  de- 
scribed by  a time-dependent  one-electron  density  operator.  In  both  cases,  the  SCF 
coupling  between  primary  and  secondary  DFs  leads  to  a primary  DF  dissipative 


Hamllionian  that  depends  upon  die  inlegrated  time  correlation  function  of  the 


vibrational  displacements  or  the  electronic  density  fluctuations. 

The  primary  DF  lime  evolution,  governed  by  the  corresponding  time  depen- 
dent  Schrddingcr  equation,  is  described  in  Chapter  4,  with  adiabatic  and  diabatic 
representation  for  the  primary  DF  electronic  basis.  The  treatment  is  then  gener- 
alized to  include  interaction  with  light,  defining  the  electron-field  representation 
[Lin  and  George  1980).  The  primary  time  evolution  operator  is  derived  by  using 
a modified  split  operator  scheme  (Feit.  Fleck  and  Steiger  1982]  for  the  dissipative 
Hamiltonian. 

Chapter  5 presents  a scmiempirical  parametrization  of  the  primary  DF  com- 
plex potential  energy  surfaces  (PESs)  and  the  transition  dipole  moment.  The 
parametrization  was  obtained  from  the  INDO  electronic  structure  calculations 
[Ridley  and  Zemcr  1973.  Zcmer  et  al.  1980]  for  a CONi  molecule  and  CO  ad- 
sorbed on  a Niio  cluster.  The  INDO  Harlree-Fock  calculations  gave  as  a result  the 
ground  state  PES.  which  was  then  rescaled  to  fit  the  previously  known  theoretical 
(Drakova  and  Doyen  1990]  and  experimental  [Hcskctt  et  al.  1985]  results.  The 
INDO  Cl  calculations  gave  the  excited  states  and  the  corresponding  transition 
dipole  moments;  the  transition  with  the  largest  dipole  moment  gave  the  excited 
state  PES  for  the  Ni(OOI)(ads)CO  system. 

Chapter  6 is  a brief  analysis  of  the  computational  aspects  of  this  work.  It  is 
primarily  concerned  with  the  numerical  time  evolution  of  wavepackcts  [Feit.  Fleck 
and  Steiger  1979,  Leforestier  et  al.  1991];  it  gives  an  outline  of  the  computaUonal 
algorithm  used  to  generate  results  presented  in  this  work,  and  a comparison  with 


other  methods.  The  flow  diagram  of  the  Fortran  77  program  used  can  be  found 
in  Appendix  B.  This  chapter  also  contains  a detailed  discussion  of  the  initial 
conditions  for  the  propagation  of  both  the  primary  nuclear  wavcfunction  and  the 
transition  dipole  amplitude,  with  the  physical  implications  of  the  choices  made. 

All  numerical  results  of  this  study  are  presented  in  Chapter  7.  they  include 
properties  such  as  expectation  values  of  the  primary  DFs  over  time,  probabilities 
of  finding  the  system  in  either  electronic  state,  the  TCFs  (autocouelation  and 
dipole  TCF),  the  desorption  lifetimes  and  lineshapes  for  several  varied  parameters 
in  calculations.  This  work's  results  are  compared  with  the  previous  ones  and 
discussed  in  connection  with  the  physical  properties  of  the  system  studied,  leading 
to  conclusions  summarized  in  the  final  Chapter  8. 


CHAPTER  2 

A PHYSICAL  MODEL  FOR  PHOTODESORPTION  OF  CARBON 
MONOXIDE  FROM  TRANSITION  METAL  SURFACES 


Qualitative  Description  of  ihe  Ni(001)(ad5)CQ  Sysiem 

The  electronic  structure  of  CO  adsorbed  on  the  transition  metal  surfaces  and 
in  particular  the  Ni(001)  surface  has  been  recently  studied  by  several  different 
approaches,  such  as  density  of  states  studies  [Sung  and  Hoffmann  1985,  Albeit 
and  Yates  1987)  or  studies  of  electronic  charge  transfer  between  the  CO  molec- 
ular orbitals  and  the  Ni  conduction  band  [Hoffmann  1988,  Zangwill  1988):  with 
experimental  results  coming  from  c.g.  electronic  energy  loss  spectroscopy  [Aki- 
moto  et  al.  19791.  Table  2.1  contains  some  relevant  data  on  the  Ni(OOI)(ads)CO 
system,  shown  in  Fig.  2.1. 

All  this  work  indicates  the  adsorption  mechanism  for  CO,  proposed  earlier 
[Hermann  and  Bagus  1977,  Sung  and  Hoffmann  1985):  its  highest  occupied 
molecular  orbital  (HOMO)  5 a donates  a fraction  of  its  population  into  the 
electronic  bands  of  nickel,  and  then  the  surface  donates  back  a smaller  fraction 
of  its  electronic  population  into  the  CO's  lowest  unoccupied  molecular  orbital 
(LUMO)  2x  — thus  a chemisorption  bond  is  formed.  As  a consequence  of  the 
interaction  with  a quasicontinuum  of  Ni  3d-4s  band  states,  the  CO  molecular 
orbitals  split  and  broaden  into  two  narrow  bands  each,  such  that  the  two  bands 
coming  from  the  5 a CO  orbital  couple  with  the  Ni  4p  upper  unoccupied  band, 
while  the  two  bands  coming  from  the  2jr  orbital  couple  with  the  3d-4s  conduction 
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band.  This  suggests  that  the  basis  set  of  electronic  states  in  which  the  system 
is  modeled  must  consist  of  modified  MOs  of  CO  — linear  combinations  of  the 
CO  HOMO  and  LUMO.  and  the  Ni  conduction  and  upper  empty  band  states.  A 
schematic  electronic  energy  diagram  for  the  Ni(OOI)(ads)CO.  showing  the  splitting 
and  broadening  of  the  CO  molecular  orbitals  due  to  the  interaction  with  the  Ni 
surface,  is  given  by  Fig.  2.2.  To  get  a further  insight  of  the  system,  Fig.  2.3 
presents  a vibrational  energy  diagram  [Parra  and  Micha  1986,  Mantca  et  al.  1990] 
for  the  Ni(001)(ads)CO,  illustrating  the  coupling  of  different  oscillating  modes  for 
the  adsorbate  with  the  phonon  bands  of  the  substrate.  Frequencies  and  symmetty 
characters  of  those  modes  are  listed  in  table  2.H. 

The  density  of  states1  crystal  orbital  overlap  population  (COOP)  results  [Sung 
and  Hoffmann  1985)  show  that  the  interaction  between  the  carbon  atom  of  CO 
and  the  adsorption  site  Ni  atom  is  bonding  in  nature  for  the  5o  orbital,  while 
it  is  antibonding  for  the  2-  orbital.  Those  results  suggest  that  the  desorption  of 
CO  is  induced  by  a transition  from  an  occupied  electronic  band  of  the  adsorbate- 
substrate  complex  that  is  a combination  of  Ni  3d(xJ  - y!)  and  Ni  3d(xy)  bands, 
into  an  empty  band  that  is  a mixture  of  Ni  4p(z)  band  and  the  carbon  2s  atomic 
orbital  which  is  antibonding  between  Ni  and  C. 


Figure  2.1: 


Sketch  of  the  Ni(OOI)(ads)CO  system 


Tabic  2J:  Experimental  data  lable  for  ihc  Ni(OOIMads)CO  system 


Properly 

Magnitude 

Reference 

Ni  conduction  band  (3d-4s)  width 

4eV 

Albert  and  Yates  '87 

Ni  Fermi  energy  level 

• 8,5  eV 

Albert  and  Yates  '87 

Ni  valence  (3p)  to  conduction 
(3d-4s)  band  gap 

> 45  eV 

Brosset  al.  '71 

Ni  conduction  (3d-4s)  to  upper 
empty  (4p)  band  gap 

< 1 eV 

Nagel  '73 

Ni  upper  empty  band  (4p)  width 

7 eV 

Nagel  '73 

Ni  Debye  frequency 

260  cm’1 

dc  Launay  'S3 

CO  binding  energy 

"•’eV 

Hen-berg  '60 

CO  equilibrium  bond  distance 

2.132  bohr 

Herzbcrg  *60 

C-Ni  binding  energy 

1.3  eV 

Heskett  et  al.  '85 

C-Ni  equilibrium  distance 

3.4  bohr 

Andcrsson  and  Pcndry 
79 

CO  (5 a to  2ir)  excitation  energy 

8.6  eV 

Akimoto  et  al.  '79 
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Figure  2.2:  Qualitative  electronic  energy  diagram  for  Ni(OOI)(ads)CO  t 


aV2itc  (cm-1) 


diagram  for  Ni(COi  )(ads 


Since  (he  Ni(OOI)(ads)CO  is  an  extended  system  with  a latgc  number  of 
degrees  of  freedom  in  the  surface,  it  is  virtually  impossible  to  study  exactly 
its  quantum  dynamics.  However,  numerous  studies  in  which  the  Ni  surface  is 
substituted  by  a finite  duster  of  Ni  atoms  (Parra  and  Micha  1986,  Drakova  and 
Doyen  1990  with  refs,  therein)  suggest  that  the  interaction  of  CO  with  the  surface 
is  quite  localized  around  the  adsorption  site. 

Due  to  the  localized  interaction,  it  should  be  possible  to  treat  the  region  around 
the  adsorbate  molecule  separately  from  the  rest  of  the  surface,  and  eventually 
apply  the  exact  quantum  dynamics  only  to  the  corresponding  degrees  of  freedom 


electronic  density  fluctuations.  The  time  dependent  Schrddinget  equation  (TDSE) 


= *]*(«))  (2.2) 
where  9(t)  is  the  total  wavefunction  of  the  system  must  be  solved  for  an  initial 
value  at  time  r = 0.  In  an  extended  system,  thermodynamical  constraints  and 
initial  conditions  lead  to  a distribution  of  the  set  of  states  {tp®}  characterized 
by  the  statistical  weights  . The  mean  value  of  a property  A at  time  r is 
therefore  given  by 

«*(*)»  = E (23) 

Two  properties  of  particular  relevance  are  the  probability  of  survival  of  state  /r 
at  time  t,  following  excitation  by  a light  pulse,  obtained  from  A = |*„(f)){*,,(t)| 


W = 


(W.K'MOI'MO) 


(2.4) 


and  the  lineshape  for  photodesorption  by  a steady  light  field,  obtained  in  the  dipole 
approximation  from  the  Fourier  transform  of 


Ul)  = (»;iP(0)tP(<)l'M0) 

^ "<*SI*S><*.WI*.(0)  125 

where  D is  a transition  dipole  operator  for  light  absorbed  during  desorption. 
Depending  on  the  wavelength  of  the  absorbed  light,  and  on  the  nature  of  the 
system,  this  transition  dipole  may  depend  on  the  adsorbate  or  the  substrate  DFs. 

The  states  'P/.ft)  must  be  obtained  in  an  approximation  suitable  to  localized 
phenomena  in  extended  systems.  One  such  approximation  is  the  self-consistent 


field  (SCF)  formulation  in  which  the  state  is  factored  into  primary  and-sccondary 
DFs  slates  i6£  and  t!>j.  This  approximation  is  limited  to  short  time  events  and  is 
most  accurate  when  used  to  calculate  expectation  values  of  additive  properties. 
In  this  case,  its  validity  may  be  expected  to  increase  because  statistical  properties 
are  being  calculated,  resulting  from  weighted  sums  over  states  and  therefore 
showing  less  susceptibility  to  approximations  to  single  states.  One  proceeds  by 
first  focusing  on  a single  state  for  given  quantum  numbers  n=(aj))  and  later 
averaging  over  the  states  0. 

In  the  self-consistent  field  (SCF)  formulation  of  the  problem,  it  is  assumed 
that  the  total  wavcfunction  of  the  system  can  be  written  as  a product  of  the  primary 
and  secondary  DFs  wavcfunctions  with  a phase  factor,  i.e. 

IMO)  = h£(0>W(0><»#X(01  (2.6) 

where  each  time  dependent  factor  is  a function  of  coordinates  of  the  correspond- 
ing region  only.  These  factors  and  the  time-dependent  phase  \ satisfy  differential 
equations  derived  with  the  time-dependent  variational  method  [Micha  and  Srivas- 
lava  1989).  The  TDSE  splits  into  two  SCF  equations  for  the  respective  DFs  of 
the  system;  omitting  the  state  indices. 

rtjyM)  = IWOWM)  (2.7) 

«^(0)  = ft(0W'W>  (2.8) 

The  effective  Hamiltonians  arc  given  by 

M0  = ffp  + M0 

(29) 

MO  = («’(i)l»p.lV-'(0)/»'(f)|^(i)) 
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and 

P.lf)  = H,  + h.[t) 

(2.10) 

MO  = W)I^W)/W)I^W> 

where  the  Dirac  bracket  notation  in  each  case  implies  integration  over  the  cor- 
responding DFs. 

A necessary  condition  for  the  validity  of  the  SCF  treatment  is  established  by 
the  nature  of  the  operators  that  describe  the  physical  process  of  interest ; these 
operators  must  depend  only  upon  the  primary  or  secondary  DFs.  but  not  on  both. 
The  case  where  the  operator  depends  on  the  primary  DFs  has  been  investigated 
for  photodissociation  with  a similar  formulation  [Srivastava  and  Micha  I 1991 ). 
In  this  case,  the  operators  arc  transition  dipoles  for  light  absorption  which  can 
depend  on  the  primary  or  secondary  DFs. 

The  SCF  calculations  provide  the  time-evolution  of  the  total  state.  If  this 
is  given  by  a wavepackct  at  the  initial  time,  corresponding  to  the  excitation 
conditions,  the  survival  probability  of  a chosen  stale  follows,  omitting  the  stale 
index  ;r.  from 

m = Pp(t)P‘{t) 

Pp(l)  = IWQ)|^(0)|i 

(^(0)|tW'(0))(^(t)|vW-(0)  (2.1 1) 

/*./,)  = l(r(0W(0)|* 

(tf’*(o)l«*(0))W(OhHO) 

and  the  spectral  lineshape  for  desorption  is  given,  when  the  transition  dipole 
depends  only  on  the  primary  degrees  of  freedom,  by 
£(()  = 

£|.(0  = M0)|Pf(0)lA.(Oltf*(O) 

{^(0)|^(0))(VK.(0|^(t))  (2.12) 

£.,0  = WQH^tO) 

W'(o)|tJ*(o)X^(t)|tf-(0) 


where  Dp(,)  = exp  (~iF,t/h)  6,(0)  exp  (ffyfh)  is  the  transition  dipole  mo- 


Fpl  = Pr  + Hi  + HrL  = F„  + £ (2. 1 3) 


(2.14) 


cue 
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dynamics  of  the  primary  DFs.  It  is  convenient  to  introduce  the  interaction  picture 
which  develops  from  the  primary  Hamiltonian, 

|^(())W  = e*p^fl’l,t)|^W)  (2.17) 

This  new  wavcfunction  is  then  introduced  into  (2.17),  which  becomes 

(MS) 


MW  = «p  ( MO  «P  (2. 19) 

The  equation  (2.IS)  is  then  formally  integrated  over  time,  and  after  transforming 
back  into  the  Schrodingcr  picture,  the  primary  wavcfunction  is  given  by 

1^(0)  =«P  (4«p<)l^(0))  P (-j«,(')'‘p(‘')l^(l'))dt' 

(2.20) 

The  procedure  is  continued  by  taking  the  solution  (2.20)  and  introducing 
it  back  into  the  last  term  on  the  right  hand  side  of  (2.18).  The  result  is 
an  intcgrodifferential  equation  for  the  primary  wavefunchon  reminiscent  of  the 
generalized  Langcvin  (GLE)  expression, 

= Hpl^(0)  + KpWll^fO))- 

The  first  term  on  the  right  hand  side  refers  to  motions  of  the  primary  DFs  in 
absence  of  the  secondary  DFs.  The  second  term  is  the  " energy  fluctuation"  tt 
its  form,  by  comparison,  contains  the  random  or  fluctuation  operator 


(2.21) 


MO  = 


(iPWItf'W) 


<-exp(-iHpt/n) 
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Finally,  the  last  term  is  the  "energy  dissipation"  term,  with  the  "memory  kernel" 
given  by 


-»IV'  exn[-;W„ ,2  23) 

■mt))  Pl  A >' 1 <^(t')i^(i'))  < 3) 

e that  Eq.(2.23)  has  a somewhat  different  structure  than 


primary  and  secondary  terms  via  projection  operators:  notably,  the  exponential 
time  evolution  operator  in  the  “memory  kernel"  is  that  of  the  primary,  rather  than 
the  secondary  DFs.  However,  it  is  obvious  that 

Mp(t,t’)  = Rr  (t)R<(t')  (2.24) 

and  that  evaluating  a statistical  average  of  both  sides  of  this  equation  over 
a distribution  of  the  secondary  stales  would  lead  to  the  familiar  fluctuation- 
dissipation  relation  (Berne  1971], 

So  far.  the  primary  integrodifferential  equation  is  formally  exact,  but  to  further 
give  its  form  it  is  necessary  to  specify  the  secondary  state.  In  the  particular  case 
of  surface  adsorption,  the  secondary  DFs  contain  a collection  of  quasiconlinuous 
excitations,  with  characteristic  frequencies  ranging  from  zero  to  values  comparable 
to  the  energy  differences  of  the  discrete  molecular  energy  levels  of  the  adsorbate. 
Therefore,  it  is  appropriate  to  calculate  the  average  fluctuation  and  dissipation 
terms  with  the  statistical  density  operator 

fyo  = £“8WJWWJ(0I  (2.25) 

where  the  toj  are  statistical  weights  for  the  set  of  initial  states  {^(O)).  e.g. 
Boltzmann  factors  for  the  temperatures  of  the  substrate,  or  nonequilibrium  weights 
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if  the  system  has  been  prepared  in  a nonstationary  distribution  of  slates.  The 
components  of  the  secondary  wavefunction  in  general  evolve  over  time  with 
ali  possible  combinations  of  the  mentioned  frequencies.  Consequently  one  can 
assume  that  the  secondary  DFs  rapidly  reach  temporary  equilibrium,  from  which 
they  can  respond  to  driving  forces;  this  leads  to  fluctuations  and  dissipation  around 
a temporary  equilibrium  which  can  be  characterized  by  e.g.  a time  dependent 
temperature. 

The  expectation  values 

4fl)« = qm) 

are  chaotic  and  rapidly  varying  functions  of  time.  The  weighted  sums  of  these 
expectation  values  and  of  their  products  can  be  expected  to  show  extensive 
cancellations,  of  the  type  described  by  models  of  stochastic  phenomena  [Kubo, 
Toda,  and  Hashitsume  1989,  van  Kampen  1981).  Therefore,  it  is  assumed  that: 

(i)  the  fluctuation  force  averages  out  to  zero  on  the  primary  time  scale,  i.e. 

= (0  = 0 (2.27) 

(ii)  the  statistical  average  of  the  memory  kernel  describes  both  instantaneous 
and  delayed  dissipations  leading  to  a temporary  equilibrium;  hence  one  writes, 
introducing  the  Dirac  delta  function,  that 

«*(M')»  = £x'54‘’,W4',>('')'  = fwytldft  - ,')  + Wr(t),(t  _ (') 

(2.28) 

where  is  an  operator  in  the  primary  DFs  only,  describing  instantaneous 
dissipation  at  each  lime,  while  describes  a remaining  persistent  dissipation 


component  with  n time  delay  given  by  a function  q analytical  in  the  variable 
This  expression  is  sufficiently  general  to  include  the  dissipation  mechanisms  that 
are  expected.  The  two  operators  to  the  right  depend  on  time  because  temporary 
thermal  equilibrium  is  allowed  for.  They  would  be  independent  of  time  if  the 
equilibrium  were  stable. 

Under  the  above  assumptions,  Eq.  (2.23)  becomes  a time-dependent  integrod- 
iffeicntial  Schrddinger  equation  with  Hermitian  and  anti-Hermitian  potentials, 

= Iff,  - !%(0JI^W>  - jff,(0  jdt1 ,(l  - (2,29) 

The  primary  Hamiltonian  is  the  sum  of  the  kinetic  energy  and  potential  energy 
operators  in  the  primary  electronic  and  nuclear  DFs, 

ff,  = ff,  + % (2  JO) 

while  the  random  force  and  memory  term  must  be  derived  from  explicit  expres- 
sions for  the  secondary  region  Hamiltonian,  and  from  the  coupling  between  the 
regions. 

A more  detailed  treatment  of  the  secondary  dynamics  is  presented  in  the 
next  chapter,  where  the  time  dependence  of  the  secondary  states  is  discussed  for 
two  types  of  excitations:  quantized  lattice  vibrations  (phonons),  and  electronic 
excitations.  Consequently,  the  secondary  Hamiltonian  can  be  expressed  as  a sum 
of  three  terms 

h.  = h;  + h;  + h?  (23  o 

where  the  first  term  refers  to  the  lattice  vibrations,  the  second  to  the  electronic 
excitations,  and  the  third  to  their  coupling.  The  coupling  Hamiltonian  can  again 
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be  expressed  as  a sum  of  vibrational,  electronic  and  mixed  terms. 

Hp,  = + H‘r,  + H (2.32) 

This  in  turn  gives  a sum  of  SCF  potentials 

hp  = hvp  + h'p  + A"  (2.33) 

in  the  primary  effective  Hamiltonian. 

To  specify  the  Hamiltonian  in  more  detail,  one  needs  to  choose  the  coordinates 
of  substrate  and  adsorbate  atoms.  Two  different  choices  are  obvious;  the  first  one 
would  be  appropriate  if  the  adsorbate  were  to  remain  always  bound  and  would 
consist  of  the  vibrational  normal  coordinates  for  the  whole  system;  the  second 
choice  introduces  the  normal  coordinates  Q,,  of  the  metal  without  the  adsorbate, 
and  the  Jacobi  (center-of-mass)  coordinates  X for  the  adsorbate  relative  to  the 
metal  surface.  In  the  present  study  of  desorption,  the  second  choice  is  more 
convenient  because  it  is  appropriate  for  the  description  of  the  asymptotic  motion 
of  the  molecule. 

Time  Correlation  Functions  and  their  SCF  Treatment 
A time  correlation  function  (TCF)  is  generally  defined  as  a projection  of  some 
time  dependent  amplitude  onto  the  value  of  that  amplitude  at  some  specified  time. 
There  exists  an  abundant  theoretical  work  demonstrating  the  relation  of  the  TCFs 
with  certain  physical  properties  of  the  studied  systems,  for  example  reaction  rate 
constants  in  terms  of  cither  Dux-flux  TCFs  (Miller  1974,  Miller  et  al.  1983) 
or  collisional  TCFs  (Micha  1979,  Micha  and  Vilallonga  1985),  dissociation  line 
shapes  in  terms  of  transition  dipole  TCFs  (Micha  1981,  Srivastava  and  Micha 
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1989],  cic.  These  relations  are  based  on  the  conversion  of  the  Fermi  golden  rule 
for  the  state  to  state  reaction  rate  (Rodberg  and  Thaler  1967.  Schiff  1969] 

^-*  = ±|HT|n>|Vm-E.,)  (2.34) 

where,  symbolically,  |m)  is  the  initial  and  |n)  is  the  final  state  of  the  system,  and 
T is  the  transition  operator;  after  summing  over  the  final  states  and  averaging  over 
the  initial  distribution,  into  an  integral  over  time  of  the  collisional  TCF  (Michn 
1981,  Micha  and  Vilallonga  1985] 

R=  J CiT(t)it 

OttW  = £u,„(n|r>(0)r(t)|n)  <235) 

f(t)  = exp  (-iffot/fi)r(0)exp 

where  is  the  Hamiltonian  excluding  the  coupling  between  initial  and  final 
slates,  and  w„  are  statistical  weights.  The  Dirac  delta  function  in  (2.34)  assures 
the  energy  conservation  in  the  system. 

The  existing  theory  must  now  be  generalized  to  include  the  dissipative  terms 
of  the  non-Hermitian  primary  effective  Hamiltonian  (2.13).  Due  to  their  presence, 
the  total  energy  of  the  primary  DFs  is  no  longer  conserved,  so  the  Fermi  golden 
rule  (2.34)  is  no  longer  valid.  The  Appendix  A presents  details  of  the  derivation 
of  the  state  to  state  reaction  rate  for  a system  with  dissipation  [BekSid  and  Micha. 
to  be  publ.|,  which  becomes 


1 (£„,  - E„f  + {Wm  - wrf/i 


(2,36) 


where  Id“)  arc  right  eigenstates  of  the  non-Hermitian  Hamiltonian  H0  = 
and  (dj|  are  their  left  counterparts.  Comparison  of  (2.34)  and  (2.36)  reveals  that 
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ihc  delta  function  in  energy  is  substituted  by  a Lorcntzian  in  energy  whose  width 
depends  upon  the  dissipation  strength  parameters  IF„  and  Wm,  that  are  dofined  as 
expectation  values  of  the  dissipative  Hamiltonian  tf^  for  the  given  eigenstates  of 
#0.  As  will  be  shown  in  the  next  section,  this  generalization  permits  the  expression 
of  the  photodesotption  lineshape  as  a function  of  the  dipole-dipole  TCF. 

The  dipole-dipole  TCF  is  a special  case  of  the  collisional  TCF  [Micha  1981, 
Mlcha  and  Vilallonga  1985);  here  flo  = Fp  + UL  and  the  transition  operator  t 
in  the  first  order  approximation  of  its  Lipmann-Schwingcr  form  (Schiff  1968)  is 
T — dpi  = —E  - D.  The  dipole-dipole  TCF  is  then 

CDD(t)  = (TD(0)|TD(<))  (2.37) 

|To(())  = [-£"«-5(<)]|«(0)) 

E(t)  = exp  (-iitit/k)£( 0)  exp  (iffgl/fi)  (2.38) 

B(t)  = exp  (—iFpt/tij  3(0)  exp  {iFpt/hj 
The  quantized  electric  field  £(t ) is  that  of  the  incident  light  and  it  oscillates  with 
time  in  the  Heisenberg  picture;  it  is  assumed  to  be  constant  over  space,  which  is 
quite  realistic  for  visible  and  near  UV  wavelengths  (order  of  magnitude  of  100 
nm)  compared  to  the  size  of  the  region  of  the  adsorption  of  the  diatomic  (order 
of  magnitude  of  0.1  nm).  Also,  as  long  the  time  dependence  of  the  radiation 
field  is  purely  oscillatory,  i.e.  E(t)  = f0exp  (-fu/gt),  which  is  the  case  for  plain 
electromagnetic  waves,  it  cancels  out  in  the  TCF  expression  (2.37)  and  the  time 
dependence  of  the  dipole  amplitude  (2.38)  is  governed  by  an  effective  frequency 
operator  defined  in  terms  of  the  primary  effective  Hamiltonian  [Srivastava  and 


Micha  I 1991),  i.e. 


IToW)  = «P  - £)«/»]  |To(0»  (2.39) 

where  E is  ihe  total  initial  energy  of  the  primary  region  of  the  system. 

The  SCF  formulation  has  been  introduced  for  the  transition  dipole  amplitude 
it  has  been  assumed  that  the  transition  dipole  operator  depends  only  on  the 
primary  degrees  of  freedom 

D = b{R,  r)  (2.40) 

so  its  lime  dependence  is  governed  only  by  the  primary  DF  Hamiltonian. 

The  dipole  amplitude  is  factored  out  [Srivaslava  and  Micha  I 1991) 

lT(0)  = KW)  ■ MO)  (2.41) 

where  the  primary  dipole  amplitude  is  given  by 

|«t))  = e*p(-iF,t)[iV(0)>]  (2.42) 

and  MO)  = |lf‘(4)>  is  the  secondary  wavefunction.  The  dipole-dipole  TCF  in 
the  SCF  formulation  then  becomes 

C0DW  = C£D(1)CJ,D(«) 

c'lt)  = «(°)l«0> 

l«(0)lfl0)){f(t)|«<f))|,/»  (2.43) 

C-U ) - (9(0)14(0) 

l(9(0)|-|(0)){,(r)|,(t))|i/2 

Another  quantity  of  interest  is  the  autocorrelation  function,  defined  as 


C(0  = (*(0)|*(f)) 


(2.44) 


Again,  the  SCF-type  factorization  can  be  introduced  and  each  factor  is  obtained 
from  the  corresponding  region  wavefunction  itself;  hence, 


where  the  separate  normalization  is  necessary  for  each  factor,  because  they  evolve 
over  time  with  non-Hermitian  complex  Hamiltonians. 

Photodcsorption  Lifetime  and  Lineshapc  in  terms  of  the  TCFs 
At  the  beginning  of  the  previous  section,  it  was  stated  that  certain  physical 
properties  of  the  system  of  interest  can  be  expressed  in  terms  of  time  correlation 
functions.  In  this  section,  it  will  be  shown  how  the  photodcsorption  lifetime  of  the 
Ni(001)(ads)CO  system  can  be  deduced  from  me  autocorrelation  function;  while 
a generalized  photodcsorption  lineshapc,  including  the  energy  dissipation  in  the 
substrate,  can  be  written  as  a temporal  Fourier  transform  of  the  dipole-dipole  TCF. 

The  absolute  value  square  of  the  autocorrelation  function  (2.59)  is  also 
known  as  Ihe  survival  probability  amplitude  (SPA).  Its  decay  rate  is  a measure 
of  the  lifetime  of  the  specified  initial  state  of  the  system.  In  this  case,  the 
decay  rate  is  computed  for  the  primary  SPA;  assuming  its  exponential  decay 
Cp(t)  = Cp(0)exp(— f/r),  it  is  defined  by  the  following  expression 


C(l)  = C'(f)C*(f) 


cm  - (r(»)ir(t» 

1 ' IW)l^(o))W)l^(0)ll/J 

c.w  _ 


(^(0)|»P(t)) 


(2-45) 


C‘U\  = (y(u)itf(t)) 
1 


(2.46) 


For  physical  systems  coupled  to  an  environment,  the  assumption  of  exponential 
decay  is  applicable  only  in  the  absence  of  such  coupling,  i.e.  at  a late  time  limit. 


Therefore,  afier  raking  a natural  logarithm  of  the  exponential  decay  law 


ln(C'(()]  = ln[C'(0)]-t/r  (2.47) 


it  becomes  obvious  that  the  lifetime  of  the  desorbing  state  can  be  estimated  as 
the  inverse  negative  slope  of  In  (Cp(i)]  in  the  asymptotic  time  limit. 

The  photodcsorption  lineshape  for  a nondissipativc  desorbing  state  was  de- 
fined as  a time  Fourier  transform  of  the  dipole-dipole  TCF  (Srivastava  and  Micha 
t 1991] 

£(“)  = 4^  |/  Cd°W  exP  M)*|  (2-48) 

where  CodIO  was  defined  by  the  eqn.  (2.37).  An  analogous  expression  must  now 
be  derived  to  include  the  dissipation  in  the  substrate. 


te  constant  (2.36).  > 


Lorcntzian  in  energy  expressed  as  a time  integral 

_{W„  - IV.) = o 

/i  J 


(E„-E„f  + (Wm-W„f/4  = 


Recalling  that  ) are  eigenstates  of  the  operator  Fo  = {Hp  - iWp/2)  + He, 
they  can  be  written  as  direct  products  of  the  eigenstates  of  the  adsorbate-substrate 
complex  and  the  eigenstates  of  the  radiation  field  Hamiltonian.  This  n 


(2.51) 
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wilh  an  analogous  expression  for  the  left  eigenvector  ($J|;  here  |u?„  ) arc  the 
eigenstales  of  the  primary  dissipative  Hamiltonian  without  the  radiation  field  with 
quantum  numbers  nr,,  and  |m*)  arc  the  photon  number  eigenstales  of  the  radiation 
field  Hamiltonian.  From  cqns.  (2.13)  and  (2.14),  it  follows  that  the  matrix 
elements  of  the  system-radiation  coupling  fipi  in  the  basis  of  photon  number 
eigenstates  are  given  by 


Introducing  (2.52)  and  (2.53)  into  (2.51),  and  defining  the  photodcsorption 
lineshape  for  the  one-photon  absorption  processes  as  a total  rate  for  the  molecular 


A!'°W  = (2.55) 

and  the  photon  number  indices  have  been  omitted.  In  analogy  with  the  analysis 
of  Appendix  A,  the  left  eigenstate  |$!,)  can  be  written  in  terms  of  its  right 
counterpart,  and  again,  if  the  dissipation  is  relatively  weak,  it  can  be  approximated 
by  W,)  in  the  eqn.  (2.55),  The  summation  over  n„  then  yields  an  identity,  and 


= YL S X)  - mu.,-0) 


(2.53) 


the  lineshape  becomes 


(2.54) 
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after  averaging  over  the  initial  distribution  of  the  adsorbate-substrate  complex 
eigenstates,  the  photodcsorption  lineshape  becomes 

L{u)  = J CDD(()ejtp(— fcrtjdl  (2.56) 

having  omitted  the  photon  indices  n*  and  m*;  with 

Cco(0  = |r«6,(0)A,(0>>  (2.57) 


CHAPTER  3 

DYNAMICS  OF  THE  SECONDARY  DEGREES  OF  FREEDOM 


Classification  of  Secondary  DF  Excitations 

As  the  secondary  region  in  this  study  is  a metal  surface,  there  are  two 
common  ways  of  excitations  that  lead  to  energy  dissipation  through  density  and 
charge  fluctuations  in  it:  lattice  vibrations  — phonons,  and  electronic  excitations. 
Consequently,  the  secondary  degree  of  freedom  (DF)  Hamiltonian  is  expressed 
as  a sum  of  two  terms 

H,  = H?  + H;  (3.1) 

and  the  coupling  to  the  primary  DFs  is  a sum  of  two  terms  as  well,  not  accounting 
for  the  electron-phonon  interactions  for  simplicity,  so  that 

h„.  = h;,+h;,  (3.2) 

To  specify  the  Hamiltonian  in  more  detail,  the  coordinates  of  substrate  and 
adsorbate  atoms  must  be  chosen.  Two  different  choices  are  obvious:  the  first  one 
would  be  appropriate  for  the  bound  adsorbate  and  would  consist  of  the  vibrational 
normal  coordinates  for  the  whole  system;  the  second  choice  introduces  the  normal 
coordinates  Q,,  of  the  metal  without  the  adsorbate,  and  the  Jacobi  (center-of-mass) 
coordinates  X for  the  adsorbate  relative  to  the  metal  surface.  In  the  present  study 
of  desorption,  the  second  choice  is  more  convenient. 


(3.3) 


%-Ea«fll*5liWi||  0.0 


wo = h : + (n‘)\HvP,m‘))/mt)m‘)) 


+E 


ig  ft,’'  and  Dp,’,  ii  can  be  written  os 

«» -?*(«•♦  J ')+ 

[i.M  (ii + i.)  + s ^(Op! + i.)  f!  t i»)  j 


mo = E Ewwiff/j  ww>/  Ewwww>  p.ii) 

^.*(0  = E Ewwi^w  i«W)/  Ewww(o)  <3. 12) 

jiv-w + *.,wi[i:w + i.w] +»,«<  = 0 (3.13) 
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and  the  Hermitian  adjoint  equation  for  tp.  Obviously,  6^'  and  blt  arc  coupled  in 
these  equations,  and  decoupling  them  is  not  simple  with  lime  dependent  coupling 
parameters  (3.11)  and  (3.12). 

It  is  useful  at  this  point  to  introduce  new  operators  for  mass  weighted 
secondary  coordinates  and  conjugate  momenta,  defined  as  linear  combinations 
of  the  creation  and  annihilation  operators 

<3'  = 7fc(iJ+s')  (3I4> 

A"Vr(^-*r)  0.'5) 

These  coordinates  and  momenta  represent  the  substrate  normal  modes  without  the 
presence  of  the  adsorbate,  and  they  satisfy  the  commutation  relation 

[£„&.]  =«V/  (3.16) 

The  secondary  DF  effective  Hamiltonian  can  be  rewritten  in  terms  of  and 
/V,  so  it  becomes 

MO  = E + y<X  + MO  <5-  + E jM(OOvf?»J  (3.17) 

The  solution  of  the  time-dependent  Schroedinger  equation  for  this  Hamiltonian 
is  a superposition  of  vibrational  stales  describing  a collection  of  normal  modes  t < 
with  occupation  numbers  n„  so  that 

wo) = Emm, oik})  au) 

Now,  instead  of  solving  the  Heisenberg  equations  for  the  creation  and  anni- 
hilation operators,  it  is  more  practical  to  solve  the  analogous  equations  for  the 


I&M  = WO 


jtPM  + «lQM  + E i [MO + MO]  QM + MO/  = o O.20) 


“ « 0f  >“ical*  «“P“  *—*  oscillator  equations  in 

£<?,(0  +"l$M  + EilMO  + M0]&(0  + wo/ = 0 0.20 


equations 


[£+nH{£!o}-°  “» 

where  by  definition  = Juj  + and  the  initial  conditions  are  chosen 

/„(0)=0;  [if Jit)  o = n„(0) 

(3.24) 

s»(0)  = i ; [igjit) o=0 

These  initial  conditions  reflect  the  fact  that  the  adsorbate  is  initially  bound  to  the 
substrate,  so  that  the  time  derivative^/,, /di)0  is  equal  to  the  initial  value  of  the 
modified  frequency  of  the  /ith  substrate  normal  mode. 

The  solution  for  the  /ith  normal  mode  (Yamashita  and  Miller  1985,  Srivastava 
and  Micha  I 1991)  is 

= <?,.(0)s,,(t)  + ^M/„(l)  - j ;f  M‘U[(<  - «')]*'  P-25) 

If  the  bilinear  coupling  terms  are  neglected  altogether,  the  time  dependence  in 
the  frequencies  disappears  and  the  functions  /,,  and  g,,  become  Ji/tfw,,/)  and 
cos(Ufit)  respectively. 

As  stated  in  the  beginning  of  this  section,  the  creation  and  annihilation  oper- 
ators £,1*  and  can  fully  describe  any  boson-type  excitations  in  the  substrate, 
such  as  phonons.  It  is  worth  noting  that  this  approach  may  be  generalized  also 
to  the  electronic  excitations  if  they  can  be  understood  as  electron-hole  pair  ex- 
citations (March  et  al.  1967);  however,  the  physical  interpretation  of  terms  in 
the  effective  Hamiltonian  is  then  very  different  because  the  lowest  term.  i.e.  ex- 
cluding couplings  to  the  primary  DFs  and  to  other  secondary  DFs'  excitations,  is 
then  linear  in  the  electron-hole  operators,  rather  than  bilinear  in  the  boson-type 


[*r4]t  = M*.]+=< 


'•»(<)  - E!«Mi*Si«»/p;wi«(i)) 


ffw- 


(3-34) 


m Va  = Vuf  m“d=  “P  of  N lime  dc 


be  defined  as  [March  cl  al.  J967J 


£ £ 


4«« -[«»««)  (3.36) 


tfie  anticommulation  rules  (3.28)  and  (3.29),  and  collecting  lerms  yields  a set  of 
coupled  time  dependent  differential  equations  for  the  one-electron  density  matrix 
elements 

= (£f  - + £ [At,f  w - Af..,p(‘)®{.p,(o] 

(3-37) 

for  each  pair  of  values  (£,  fc1). 

To  the  lowest  order,  neglecting  the  coupling  altogether,  the  above  set  of 
equations  has  a simple  oscillatory  solution 

= «P  (-i(cf-£c,)tje™(0)  (3.38) 

When  the  coupling  is  included,  the  equations  of  motion  for  the  one-electron 
density  matrix  elements  can  be  solved  using  the  metal  dielectric  function  if  the 
coupling  is  weak  [Kittel  1963),  or  e.g.  by  temporal  linearization  [Micha  and 
Runge  1992)  if  the  coupling  is  strong.  The  operator  solution  is  given  in  terms  of 
the  time-evolution  operator  for  this  potential  as 

M‘)  = V,(t,0)p,(0)U,(t,0)'  (3.39) 

As  in  the  case  of  phonons,  initial  conditions  must  be  specified  taking  into  account 
that  the  adsorbate  is  bound  to  the  substrate  at  the  initial  lime. 

Form  of  the  Primary-Secondary  DF  Coupling  Hamiltonian 
In  this  section,  the  primary-secondary  DF  coupling  Hamiltonian  is  derived 
in  lerms  of  the  solutions  to  the  secondary  DF  dynamics  for  both  phonons  and 
electrons.  Then,  using  the  definitions  from  Chapter  2,  the  form  of  the  primary 


K = E m®)  E [i;jX)q. + W)<8]  mi  cm<q 


(<?4<))fl  = «0(OI&WJ(O)/WJ(OWJW)  = <£{t)  (3.4i) 


hm,  - <')*'  (343) 

"{^yifcn  - j M,.(,  _0i,J 


Finally,  the  expectation  value  of  the  coupling  Hamiltonian  (3.7)  can  be  written 
in  terms  of  (3.41)  and  (3.42).  so  it  becomes 

(3.45) 

hh(X,<)  = Y. 


the  metal  lattice  vibrations  will  be  oscillating  rapidly  due  to  the  summation  over 
normal  mode  contributions. 

The  matrix  elements  of  couplings  between  different  electronic  states  of  the 
adsorbate  come  from  the  momentum  couplings  of  the  Bom-Oppcnheimer  theory, 
after  the  transformation  from  the  adiabatic  to  the  diabatic  representation;  they 
are  likely  to  be  small,  except  at  crossings  or  avoided  crossings,  and  arc  likely 
to  depend  only  weakly  on  lattice  displacements;  therefore  a good  assumption 
is  that  the  coupling  matrix  is  diagonal  in  the  electronic  labels  and  linear  in  the 
displacements,  giving 

= E W(*»4?(*.t)W(*)l  (3.46) 

*//(*.<)  = (3.47) 

where  the  sccondaty  slate  index  0 has  been  restored.  The  vibrational  excitations 
of  the  substrate  would  therefore  contribute  to  the  random  forces  and  dissipation  in 
the  electronically  adiabatic  motion  of  the  adsorbate,  but  would  not  contribute  to 
electronic  transitions.  Expressions  similar  to  the  one  above  have  been  frequently 
used  in  the  literature  to  model  coupling  to  a vibrating  lattice  [Messina  and  Coalson 


*:-Ei  Mtmlt.  «)W(*)I 
%(*,<)  -EE»ttWs3WE»MW 
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leads  10  a very  rapidly  varying  iracc.  so  lhat  ihc  fluciualion  and  raeraoiy  lerms 
(2.27)  and  (2.28)  of  Ihe  primary  effective  Hamillonian  should  show  the  lime 
dependence  typical  of  Brownian  motions. 

The  treatment  above  can  be  extended  to  include  interband  transitions.  In 
this  case  the  quasipatticle  states  must  be  labelled  by  the  band  indices  b = C.U. 
for  the  conduction  and  unoccupied  bands  respectively,  in  addition  to  the  quasi- 
momentum label.  A new  index  \ = {b,kj  can  be  used  to  write  the  more  general 
matrix  elements  of  the  coupling  and  the  density  matrix.  The  double  sums  in 
h)j  can  be  rewritten  as  the  sum  of  terms  relating  to  the  conduction  intraband 
transitions,  plus  terms  relating  to  interband  C — U and  C—  V transitions.  The 
SCf  coupling  in  the  primary  states  can  be  written,  restoring  the  8 index,  as 

W = E W(*)|  (3-51) 

ru'W-iil/EAw 

in  terms  of  a normalized  electron  density  matrix  7. 

The  potential  matrix  elements  Ufa  derive  from  both  the  diabatic  coupling 
between  electronic  and  nuclear  motions,  and  from  the  Coulomb  interaction  among 
electrons.  They  are  larger  when  the  adsorbate  electronic  state  is  the  same  / = J, 
because  then  there  is  a larger  overlap  of  the  total  electronic  states.  For  / # J 
there  is  also  an  interaction,  with  a longer  range  than  vibrational  couplings  because 
ihc  present  electronic  couplings  have  a Coulombic  origin. 

So  far  the  treatment  of  the  secondary  DFs  has  considered  the  boson-like 
excitation  operators  (fly}  = (Qy.Py)  for  the  lattice  vibrations  and  {fly}  = 
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{sjfi /}  for  the  electron-hole  excitations.  To  complete  the  treatment  of  electronic 
excitations  in  the  substrate,  one  should  also  consider  the  coupling  of  e-h  excitations 
with  lattice  vibrations.  This  can  be  done  introducing  collective  excitation  operators 
{Bj}  = (“lai'Qji  , . which  combine  electron-hole  and  phonon 

operators,  and  proceeding  as  before.  However,  this  shall  not  be  done  for  sake 
of  simplicity;  instead,  the  effective  primary  Hamiltonian  will  be  constructed  for 
the  previous  two  cases. 


Form  of  the  Primary  DF  Energy  Dissipation  Operator 


From  the  definition  in  Chapter  2 for  the  energy  fluctuation  and  memory  terms, 
introducing  now  the  state  index  0 for  the  secondary  states 


WONgffl)  e*p(~l/W,‘)  (3’53) 


Jii*V)  = 


WmW)  ' 


mnw)) 


(3.54) 

These  operators  may  be  obtained  for  substrate  phonons  or  for  e-h  excitations  and 
can  be  averaged  over  Ihe  statistical  distributions  of  secondary  states  0.  In  both 
cases  they  arc  given  in  terms  of  the  expectation  values  of  the  collective  motion 
operators,  by 


A?1  = (*#(') 

*“(<)  = (s>(0)fl 


(3.55) 
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The  statistical  averages  of  fluctuation  and  dissipalion  terms  are  then 
*pW  = £«W  t)))Hi(X)exp(-iHptlh) 

km -E-rim 

*«•(</)  = - t')/h\H,(X) 

«6j(0V(I')»  = E’“^(‘)^(O 

These  expressions  are  similar  to  Ihe  ones  lhai  appear  in  Ihc  litcraiure  [Kubo 
et  al.  1989,  Blum  1981]  and  can  be  analyzed  analogously.  The  fluctuation  and 
dissipalion  terms  of  the  effective  Hamiltonian  include  the  average  of  the  collective 
motion  operators  and  the  time-correlation  of  pairs  of  them,  respectively.  The 
assumption  of  Brownian  motion  for  the  adsorbate  interacting  with  the  substrate 
means  that  the  substrate  DFs  oscillate  rapidly  and  irregularly  on  the  lime-scale  of 
Ihe  adsorbate  motions,  so  that  the  substrate  rapidly  reaches  temporal  equilibrium. 
Therefore,  at  each  time  one  can  write 

«*/(<)».  = 0 (3.58) 

and  that  the  time-correlation  function  (TCF)  is  independent  of  the  lime  r, 

«‘j(')WO».  = «M0)V(r))>,  (3.59) 

where  t = t'  -t  and  the  subindex  r indicates  that  the  average  is  over  a temporary 
equilibrium  distribution. 

This  time-correlation  function  may  be  analyzed  into  its  Fourier  components  to 
discover  how  the  secondary  DF  frequency  range  will  shape  the  lime-correlation. 


"w*‘> = «■**  E^(0;0E  W)^v(^) 


after  long  limes.  This  makes  it  possible  lo  deal  wilh  substrates  described  by  a 
time-dependent  temperature,  in  accordance  with  the  recent  experimental  results 
(Cavanagh  et  al.  1993,  Gcrmer  et  al.  1993,  Prybyla  ct  al.  1990  and  1992). 
More  detailed  expressions  for  the  dissipation  potential  can  be  found  from  the 
derivations  for  phonon  and  c-h  excitations. 

For  phonons,  from  the  corresponding  section,  using  the  linear  SCF  potential, 
= E W(*)>‘»(*.  tmm\  (3.64) 

Afj(*.*)  = Eff"(-*)rf(0  (3.65) 

^=Ylwh{x^m[x)m{x)  i 

W!AX,  t)  = 4rh-‘  £[£  HlK(X)H?J(jC)\  j dr  {(,„(<), „(r))),  ^ 


of  the  vibrational  displacement  averages  q„(l)  . These  can  be  replaced  [cqn. 
(3.25))  io  obtain 

/ * J dr  g.{r) 

cc  (3.67) 

»««« 1 j dr  /,(t) 

From  the  section  on  phonon  dynamics,  it  is  also  apparent  that  the  average  over 
all  times  appearing  in  the  dissipation  term  would  give  zero  if  the  effect  of  the 
adsorbate  on  the  secondary  region  were  neglected,  because  f„(l)  would  then  be 
combinations  of  sine  and  cosine  functions  and  their  integrals  would  vanish  except 
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for  a vanishing  lattice  frequency.  The  only  terms  in  qr(t)  which  contribute  to 
the  instantaneous  dissipation  come  from  the  self-consistent  coupling  that  has  been 
introduced,  which  alters  the  oscillations  in  the  secondary  DFs.  Dissipation  through 
lattice  vibrations  would  be  absent  were  it  not  for  the  self-consistent  treatment. 
For  electron-hole  excitations,  using  results  from  the  corresponding  section. 

= E l^l)hff(x,  <)«$(*)!  (3.68) 

0-69) 

This  gives,  in  analogy  with  the  phonon  case 

K (3.70) 

Cu..x»i»(0;l)=  J d-r 

Hence,  the  dissipation  potential  depends  on  the  integrated  TCF  of  electron  density 
fluctuations,  which  con  relate  to  intraband  or  interband  transitions.  For  both 
intraband  and  interband  dissipation,  the  TCFs  appear  integrated  over  all  times 
rn  the  instantaneous  dissipation  term.  This  term  would  vanish  identically  if  the 
effect  of  the  adsorbate  on  the  secondary  DFs  were  ignored,  because  the  density 
fluctuations  would  then  oscillate  harmonically.  The  scir-consistcnt  field  treatment 
of  this  work  gives  a nonzero  result  because  it  includes  in  the  secondary  effective 
Hammoman  a term  describing  electron  state  couplings  due  to  the  presence  of  the 


adsorbate. 


CHAPTER  4 

TIME  EVOLUTION  OF  THE  PRIMARY  DEGREES  OF  FREEDOM 

TD  Schrodingcr  Equation  for  the  Primary  DFs 


This  chapter  studios  the  primary  DFs  quantum  dynamics  of  the 
Ni(100)(ads)CO  system.  It  is  governed  by  the  time  dependent  Schrfldingcr 
equation  (TDSE)  with  the  primary  DF  effective  Hamiltonian,  derived  from  the 
SCF  treatment  of  Chapter  2, 

iJl^W>  = W(0)  (4.1) 

where  ft  = 1 here  and  henceforth.  The  wavcfunction  tPU)  is  a function  of  all  the 
primary  coordinates  and  time.  It  fully  describes  the  dynamics  of  the  adsorbate 
in  the  presence  of  the  substrate,  which  is  modeled  as  a dissipative  medium.  The 
primary  effective  Hamiltonian,  from  Chapter  2,  is 

tut 

where  R is  the  distance  between  the  CO  center  of  mass  and  the  surface,  and  r is 
the  C-0  distance.  Due  to  the  mass  difference  between  electrons  and  atomic  nuclei, 
the  time  scales  for  electronic  motion  are  much  slower  than  those  of  the  nuclear 
motion,  so  it  is  reasonable,  in  the  spirit  of  the  Born-Oppcnhcimcr  approximation, 
to  separate  the  nuclear  kinetic  energy  (the  first  two  terms  in  brackets)  from  the 
rest  of  the  Hamiltonian  — the  complex  potential  energy.  The  equation  (4.2) 
implies  that  the  imaginary  potential  depends  only  on  the  primary  coordinates  and 
not  their  conjugate  momenta.  This  is  justified  by  the  physical  meaning  attributed 
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(o  the  imaginary  potentials  — broadening  of  ihe  discreie  quantum  energy  levels 
of  the  primary  DFs.  The  form  of  Ihe  imaginary  potential  was  derived  in  detail  in 
Chapter  3 for  the  secondary  DFs  describing  vibrational  and  electronic  excitations 
of  the  substrate,  and  the  above  assumption  is  confirmed. 

It  is  suitable  to  choose  a basis  of  electronic  states  (parametrically  dependent 
upon  nuclear  variables)  in  which  to  expand  the  primary  wavefunction;  Ihe  coef- 
ficients of  expansion  are  then  functions  of  both  nuclear  variables  and  time 

llPM)  = £W(*,r.fM(K1'0>  (4.3) 

and  the  primary  effective  Hamiltonian  is  given  by 

= E W(*. r))(*H  + Vu  - iW,il2)Wj(R,  r)|  (4.4) 

where  each  Ktj  is  a matrix  element  of  the  primary  kinetic  energy  operator  on 
nuclear  variables,  and  Vjj  -iW/j/2  is  a matrix  element  of  the  complex  potential 
energy  function  of  the  nuclear  variables.  The  electronic  states  |^5(/i,r))  for  the 
case  of  Ni(ads)CO  are  the  molecular  orbitals  (MOs)  of  the  adsorbate-substrate 
complex,  and  not  of  the  isolated  CO. 

The  expansion  (4.3)  is  then  introduced  into  the  TDSE,  so  it  can  be  written 


= [*  + V(V>  - jW(fi,r)]r(')  (4.5) 

where  K is  the  kinetic  energy  operator  on  nuclear  variables,  an  N by  N matrix. 
V and  W are  the  N by  N matrices  of  the  complex  potential  energy,  and 
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The  elements  of  the  kinetic  energy  matrix  are 

Mu  ' >55 

where  the  second  term  in  each  square  bracket  represents  the  momentum  coupling 
between  the  electronic  basis  states. 

The  elements  of  the  real  potential  energy  matrix  are  given  by 

(v(«.r)),,j  = MHV)l?Wr)W(fl,r))  (4.7) 

and  they  represent  all  Coulomb-type  interactions  between  the  electronic  basis 
states.  Finally,  the  imaginary  potential  matrix  elements  are  given  by 

(W  (B.f))W  = (<ff(*,r)|H'(K.r)|tf5(/l,r))  (4.8) 

Again,  as  in  Chapter  2,  the  above  formalism  must  be  expanded  to  include 
interaction  with  light.  Since  the  radiation  field  is  treated  quamum-mcchanically 
[Louisell  1973,  Srivastava  and  Micha  I 1991),  it  is  convenient  to  define  a new 
electronic  basis  as  a direct  product  of  the  basis  of  adsorbate-substrate  electronic 
orbitals  and  the  photon  number  wavefunctions  labelled  as  |n*),  so  that 

14  "a)  = (4,9) 

where  Itj  is  the  photon  number.  This  basis  constitutes  a so-called  electron-field 
representation  (Lin  and  George  1980).  This  study  of  the  Ni(001)(ads)CO  system 
interacting  with  light  includes  electronic  excitations  and  one-photon  processes, 
so  the  primary  basis  states  are  |j,  1*)  (ground  electronic  state  of  the  adsorbate- 
substrate  complex  with  one  photon  present)  and  |e,  0*)  (excited  electronic  stales 
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with  no  photons).  In  this  basis,  the  matrix  of  the  radiation  field  Hamiltonian  is 


(4.10) 


while  the  system-field  coupling  Hamiltonian  matrix  is 


where  D(R,r)  is  the  transition  dipole  moment. 

The  next  task  is  to  find  a suitable  basis  set  of  electronic  states  to  evaluate 
these  matrices,  i.e.  to  choose  an  electronic  representation.  There  is  an  infinite 
number  of  possibilities,  but  the  choice  is  made  to  facilitate  the  numerical  solution 
of  the  problem. 


Of  all  possible  electronic  representations,  two  are  especially  useful  — in  each 
of  them  one  of  the  terms  of  the  Hamiltonian  matrix  is  diagonal.  The  first  one 
eliminates  the  electronic  coupling  between  the  different  basis  states  through  the 
potential  energy  matrix  and  it  is  called  the  adiabatic  representation  [Smith  1969, 
Baer  1975,  Olson  and  Micha  1982,  Heather  and  Meliu  1987,  Fernandez  and 
Micha  1991,  Cohen  and  Micha  1992).  Hence 


The  other  possibility  is  the  representation  in  which  there  is  no  momentum 
coupling  between  different  electronic  basis  states,  and  is  called  the  strictly  dinbatic 
representation.  One  writes 


Electronic  Representations 


\2McoOH2  + 2/ico'i  1 


(4.13) 


where  Mco  is  the  total  mass  of  CO.  pc0  is  its  reduced  mass,  and  I is  the  identity 
matrix.  In  this  representation,  the  matrices  of  the  complex  potential  energy  are 
nondiagonal  and,  for  a particular  case  of  two  electronic  states,  their  explicit  forms 


V.j(K,r)j 

Vn{R,r)\ 


(4.14) 


w(4  _ ['’'''ll ((*.’■)  Wn{R,r)' 

K(«.r)  W-alR.r) 
The  physical  model  in  this  study  is  developed  in  the 
the  momentum  couplings,  with  a form  usually  unknown,  a 


(4.15) 


! replaced  by  functions 


the  resulting  electronic  couplings  between  different  basis  states  contribute  to  both 
the  real  and  the  imaginary  part  of  the  primary  DFs’  Hamiltonian.  These  couplings 
have  two  physical  causes:  interaction  between  the  primary  and  secondary  DFs 
charge  distributions  (predominantly  the  adsorbate  dipole  moment  with  the  image 
dipole  in  the  substrate),  and  the  momentum  couplings  between  the  primary 
electronic  states. 

The  elements  of  the  real  and  imaginary  potential  matrices  (4.14)  and  (4.15) 
arc  all  functions  of  the  primary  coordinates  only,  therefore  both  V<4>  and  \VMI  arc 
real  and  symmetric  matrices  (so  that  Va  = V2/  and  similarly  for  the  imaginary 
part).  This  means  that  they  can  be  diagonalized  by  orthogonal  transformations  as 
a special  case  of  diagonalization  of  Hermitian  matrices  by  unitary  transformations 
[Arfkcn  1985).  An  important  point  here  is  that  these  unitary  transformations 
are  different,  which  will  have  to  be  taken  into  account  later  when  the  primary 
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DF  lime  evoluiion  operator  is  consmicicd.  Also,  they  will  depend  on  (he  initial 
conditions  of  propagation,  i.e.  on  whether  Ht  and  Hpl  should  be  included;  if  so. 
the  imaginary  part  of  the  Hamiltonian  matrix  becomes  a sum  -iW(d,/2+HpL. 

For  the  real  part  of  the  potential,  the  diagonalization  of  V<dl  is  in  fact  the 
transformation  into  die  adiabatic  representation,  as  defined  for  real  potentials 
[Smith  1979,  Baer  1975,  Olson  and  Micha  1982];  it  is  given  by 

V<*>(R,r)=LyV(d,(*,r)-4  (4.16) 


where  the  transformation  matrix  Lv  is  an  element  of  the  SO(N)  algebra  [Fernandez 
and  Micha  1991]  (in  this  particular  case  N=2),  and  can  be  parametrized  with  a 
rotation  angle  that  is  a function  of  primary  coordinates, 

ZVMt) 


= ; 

Then  the  transformation  man 
Lv  = 


Lv.t(fi,r) -/„(/!,  r)J 


(4.17) 


Mow,  an  analogous  transformation  is  defined  for 
ix  if  the  interaction  with  light  is  omitted,  it  is  defined  by 
WD(rt,r)  = LW  • W<d>(/i,r)  ■ L[v 


(4.18) 

imaginary  potential 


r «*<(/!, r)  *,«*,,)  1 

[-sin  ((71,  r)  «*<(*, r)J 


(4.21) 
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It  must  be  stressed  here  that  the  diagonalization  of  the  imaginary  potential  operator 
does  not  generate  a representation  in  the  physical  sense,  as  it  does  not  eliminate 
any  physical  coupling,  but  is  merely  a convenient  mathematical  procedure  that 
will  be  used  in  the  construction  of  the  primary  DF  time  evolution  operator. 

The  physical  model  of  this  study  is  constructed  in  the  diabatic  electronic  rep- 
resentation; the  momentum  couplings  in  the  adiabatic  representation,  which  are 
large  only  near  potential  crossings  or  avoided  crossings,  are  replaced  by  smooth 
coupling  functions  of  coordinates  in  the  diabatic  representation  obtained  from 
physical  considerations;  the  resulting  electronic  couplings  between  basis  slates 
contribute  to  both  the  real  and  the  imaginary  part  of  the  coupling  Hamiltonian  be- 
tween primary  and  secondary  DFs.  When  the  couplings  are  instead  available  from 
calculations  for  surfaces  or  clusters  in  the  adiabatic  representation  [Hcad-Gordon 
and  Tully  1992),  the  diabatic  matrix  elements  follow  by  the  above  transformation. 

Time  Evolution  Operator  with  the  Split  Operator  Scheme 

The  formal  solution  to  the  TD  Schrbdingcr  equation  for  the  primary  DFs  can 
be  written  immediately 

|t^(i  + At))  = exp  |^n(t))  (4.22) 

V"(t  + At)  = exp  [-iAt^K  + V - |w^"(t)  (4.23) 

The  time  evolution  operator  is  therefore  an  exponential  of  a sum  of  two 
operator  terms  — kinetic  and  potential,  and  the  complex  potential  term  is  a sum 


of  two  terms  itself.  For  any  two  operators  A and  8.  the  exponential  of  their  sum 
can  be  expressed  approximately  as  the  Weyl's  formula 

exp  (a  + fl)  = exp  (a)  exp  (b)  exp  ^-j[a,b]  j + h.o.t.  (4.24) 
This  expiession  comes  from  the  Magnus'  expansion  [Magnus  1954]  for  the 
product  of  two  exponentials  of  operators,  truncated  at  the  first  order  correction. 
It  is  the  basic  idea  of  the  so  called  split  operator  scheme  [Feit,  Fleck  and  Steiger 
1982.  Leforestier  et  al.  1991],  where  the  time  evolution  is  executed  first  by  half  a 
step  with  the  potential  propagator,  then  by  a full  step  with  the  kinetic  propagator, 
and  finally  again  by  half  a step  with  the  potential  propagator.  As  a consequence, 
the  residual  error  is  quadratic  in  the  time  step,  rather  than  linear  as  it  would  be  by 
simply  approximating  the  exponential  of  the  sum  as  the  product  of  exponentials. 
To  account  for  the  complex  potentials  that  appear,  the  same  idea  will  be  applied 
on  the  potential  propagator  itself. 

Taking  advantage  of  the  relatively  simple  way  of  diagonalizing  the  potential 
and  kinetic  energy  matrices  for  a small  electronic  basis  (N=2).  the  solution  to  the 
TD  Schrodinger  equation  can  be  written  with  the  split  operator  scheme,  such  that 
all  three  exponents  in  the  time  evolution  operators  are  diagonal  in  the  respective 
bases  of  electronic  slates.  The  solution  for  real  potentials  reads  [Femindez  and 
Micha  1991] 

v»“(t  + Af)  = Lv*Uy^(A(/2)  • Lv  - 

t (.1  <«S) 

I'v  • U'“’(At/2)  • LV  ■ *»(i) 


"’(At/2)  = exp  (-iAlVw/2) 


(4.26) 


is  the  potential  energy  propagator  and 


U”>(A()  = exp  (-.AlK,d>)  (4.27) 

is  the  kinetic  energy  propagator. 

Starting  from  the  equation  (4.23),  and  making  use  of  the  associativity  of 
matrix  addition,  i.e.  H<d>  = (K<d>  + V<d>)  - <W«*>/2;  the  imaginary  potential 
may  be  split  first  and  then  the  real  potential  and  the  kinetic  terms  may  be  split 
inside  of  it.  The  result  for  the  primary  DFs'  wavcpnckct  propagation  in  that  case 
becomes  a generalization  of  the  time  evolution  expression  (4.25). 

W + At)  = (L^  • Uwd(A</2)  • LW|- 
[L'v  • U(v’(At/2)  ■ LV]  • U(k’(A0- 

, r , (4-28> 

[Ly  - Uy  (At/2)  • LV]' 

■p-w  ■ Uwd (At/2)  • Lw)  • v>°(l) 
where  the  imaginary  potential  propagator  is  given  by 


Uwn(A(/2)  = exp  (-AtWD/4)  (4.29) 

where  Wo  is  the  diagonalized  imaginaty  potential  matrix,  and  the  other  two 
remain  the  same  as  before. 

The  expression  (4.28)  for  the  time  evolution  over  coupled  complex  potential 
energy  surfaces  can  be  called  the  “double  split  operator  scheme",  as  the  original 
split  operator  scheme  (Feit,  Reck  and  Steiger  I982|  is  employed  on  the  whole 
propagator  to  separate  first  the  imaginary  potential  contribution,  and  then  again 
on  the  primary  Hamiltonian  propagator  to  separate  its  kinetic  and  real  potential 
contributions. 


CHAPTER  5 

POTENTIAL  ENERGY  FUNCTIONS  AND  THEIR 
COUPLING  FOR  PRIMARY  DEGREES  OF  FREEDOM 


Primary  DF  Potential  Energy  Surfaces  and  their  Couplings 
The  physical  model  for  photodesorplion  of  dialomic  molecules  from  transilion 
metal  surfaces  was  introduced  and  developed  in  the  previous  chapters.  Now  it 
must  be  specified  for  the  electronic  potential  energy  surfaces  (PESs)  as  functions 
of  primary  DFs;  they  are  later  used  for  studying  the  quantum  molecular  dynamics 
of  the  Ni(OOI)(ads)CO  system. 

The  ground  state  PES  must  describe  a bound  state  of  the  adsorbate  on  the 
substrate  surface,  therefore  it  must  have  a minimum  for  some  equilibrium  distance 
R,  between  the  adsorbate’s  center  of  mass  and  the  adsorption  site.  In  the  simplest 
case  [Menzel  and  Gomcr  1964,  Redhead  1964,  Hermann  and  Bagus  1977],  the 
whole  system  is  treated  as  a quasi-diatomic  and  the  ground  state  PES  may  be 
defined  e.g.  as  a Morse  potential  in  the  distance  R between  the  surface  and  the 
center  of  mass  of  the  adsorbate  that  is  treated  effectively  as  an  atom 


The  form  of  the  excited  electronic  state  is  expected  to  vary  for  different 
systems.  It  may  be  a bound  state,  so  it  can  also  be  described  by  a Morse  potential 
with  a constant  displacement  to  account  for  the  asymptotic  energy  difference 
t/"(R)  = U _ 2e-n.W-fl.)j  + AV- 


(5.2) 
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In  this  case,  the  desorption  is  possible  through  a DIMET  process  discussed  in 
Chapter  2,  an  example  being  the  Pt(ads)NO  system  [Gadzuk  1991,  Misewich  el 
al.  1992],  On  the  other  hand,  the  excited  electronic  state  may  be  repulsive  in 
K.  as  assumed  by  the  MGR  model,  so  one  may  write  the  MGR-type  repulsive 
potential  in  the  reaction  coordinate  R as 


£/«"7R(fl)  = Vie—'*  + AK,  (5.3) 

Such  systems,  Ni(ads)CO  being  one  of  them,  evolve  over  the  excited  state  PES 
after  a Franck-Condon  type  transition  from  the  ground  state  and  thus  directly 
desorb  from  the  substrate. 

Clearly,  the  potential  (5.1)  docs  not  include  any  internal  properties  of  the 
adsorbate.  For  diatomics,  it  is  reasonable  to  expect  that  the  internal  vibration  will 
be  excited  in  the  desorption  dynamics;  however,  the  typical  diatomic  bonds  are 
much  tighter  than  the  typical  chemisorption  bonds,  for  example,  the  desorption 
energies  for  CO  adsorbed  relative  to  the  Ni  surfaces  [Heskctl  et  al.  1985]  are 
around  1.3  eV.  compared  to  C-0  dissociation  energy  of  1 1.3  eV  (Hcrzbcrg  I960]. 
This  means  that  the  dissociation  of  the  adsorbate  is  unlikely  for  the  processes 
of  interest  in  this  work  (light  wavelength  of  approx.  250  nm  equivalent  to 
approximately  5 eV  energy  — substantially  less  then  the  dissociation  energy  of 
CO),  so  the  internal  vibration  of  the  adsorbate  can  be  described  as  a harmonic 
potential  well  in  the  diatomic  distance  r,  around  the  equilibrium  value  rt 

(U> 

A further  refinement  of  the  Mcnzcl-Gomer-Redhcad  (MGR)  model  is  achieved 
by  taking  into  account  the  interaction  of  the  permanent  dipole  moment  of  the 
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adsorbate  (such  as  CO)  with  the  induced  dipole  moment  that  its  presence  creates 
in  the  substrate.  It  is  assumed  that  the  adsorbate  permanent  dipole  moment  can 
be  expanded  around  the  equilibrium  internal  distance,  up  to  the  linear  order  in 
r,  it  then  becomes 

rf»(r)  = [‘isO  + rf,i(r-r,))nr  (5,5) 

The  aforementioned  interaction  is  then  described  as  the  interaction  of  the  above 
dipole  with  its  image  in  the  surface,  so  the  potential  energy  is 

which  is  just  one  half  of  the  interaction  of  two  antiparallel  dipoles  separated  by 
R.  The  factor  1/2  comes  in  due  to  the  fact  that  one  of  them  is  an  image  dipole 
(Jackson  1975]  and  therefore  the  energy  density  of  that  interaction  is 


where  is  the  induced  polarization  of  the  medium  and  E„  is  the  external 
electric  field. 

Finally,  the  ground  state  PES  is  obtained  by  adding  up  the  contributions 
(5.1),  (5.4)  and  (5.6),  so  the  final  result  is 

V„(rt,r)  = V* {e-7‘»,l*-«,-7<r-r,)I  _ 2e-n1|fl-fi1-,(r-r()|  jx 

x[I  + b,(r  - r,))+  (5,8) 

where  7 = mc/(mc  + mD)  is  the  mass  factor  that  appears  because  of  the 
definition  of  S (Ni  to  CO  center  of  mass  distance,  rather  than  the  Ni-C  distance). 
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The  first  term  of  this  potential  function  is  a Morse  potential  in  the  Ni-C  distance 
describing  a relatively  weak  chemisorption  bond  between  Ni  and  CO,  multiplied 


Rg  and  r,  denote  the  minimum  of  the  potential.  The  second  term  of  the  ground 
state  PES  is  a harmonic  potential  in  r describing  a strong  chemical  bond  between 
C and  0,  allowing  for  vibrations,  but  not  the  dissociation  of  the  adsorbate.  The 
third  term  describes  the  interaction  between  the  permanent  dipole  moment  of  CO. 
expanded  up  to  the  first  order  in  r,  and  its  image  dipole  induced  in  the  substrate; 
a switching  function  is  multiplied  into  it  to  provide  a cutoff  at  short  distances. 

The  ground  electronic  state  of  the  Ni(OOI)(ads)CO  system  has  been  studied 
thoroughly  (Sung  and  Hoffmann  1985,  Albert  and  Yates  1987,  Zangwill  1988); 
parameters  for  the  ground  state  PES  given  by  the  cqn.  (5.8)  are  presented  in  table 
ST,  with  references.  The  parameters  of  the  dipole  interaction  and  the  coupling 
botween  R and  r were  obtained  from  electronic  structure  calculations  — details 
are  given  in  the  final  section  of  this  chapter. 

The  imaginary  part  of  the  ground  state  potential  comes  from  the  coupling  of 
the  CO  5(7  orbital  with  the  substrate  lattice  vibrations  and  intraband  electronic 
transitions  within  the  Ni  conduction  band,  as  shown  in  Chapter  3.  This  coupling 
must  vanish  for  large  values  of  R.  as  the  adsorbate  remains  in  its  ground  slate 
asymptotically.  The  chosen  parametrization  is 

W,[n,r)  = W°exp{-2a,(R  - 7r  - «,)) . [1  + 6,(r  - r,)|  (5.9) 

i.e.  the  functional  dependence  is  the  same  as  that  of  the  repulsive  term  of  the 
Morse  potential  in  V,( R,r),  with  a separate  strength  parameter.  As  was  shown 


in  Chapter  3.  the  dissipation  strength  parameter  IP,0  is  related  to  the  density  and 
charge  fluctuations  in  the  substrate;  however,  in  this  study  it  is  left  as  a variable 
parameter  in  calculations. 


In  analogy  with  the  ground  electronic  slate,  the  dipole-dipole  term  and  the 
:rna!  vibration  harmonic  potential  are  added  to  the  exponentially  decaying 


repulsive  potential  in  R,  so  the  complete  excited  state  PES  is  given  by 
K(*,r)  = V} exp  - yr)]  • [1  + 6„(r  - re)J  + AVJ+ 

, pfif,  .2  fdeo  + dei(r  — r-)l2  f I 1 ) (5.10) 

+~f(r + 

The  first  term  is  a repulsive  exponential  function  in  the  Ni-C  distance,  describing 
the  antibonding  between  Ni  and  C.  again  with  coupling  of  the  nuclear  coordinates. 
The  second  term  is  the  asymptotic  energy  difference  between  the  Ni  conduction 
and  lowest  unoccupied  bands.  The  last  two  terms  are  analogous  to  the  corre- 
sponding terms  in  the  ground  state  PES.  The  parameters  of  the  excited  stale  PES 
(5.10)  were  obtained  by  fitting  the  results  of  the  electronic  structure  calculations, 
explained  in  a later  section. 

The  imaginary  part  of  the  excited  state  potential  is  again  parametrized  as  a 
repulsive  exponential  function  in  the  Ni-C  distance 

Wc(R,t)  = VP® exp  [-<*,(«  - 71-)]  ■ (1  + 4e(r  - re)|  (S.l  I) 
Again,  the  dissipation  strength  parameter  \V,°,  related  to  the  substrate  dynamics, 
is  left  as  a variable  parameter. 

It  is  apparent  that  this  model  of  adsorption  docs  not  take  into  account 
the  overall  rotation  of  the  adsorbate,  it  is  assumed  that  its  symmetry  axis  is 
perpendicular  to  the  surface  and  thus  no  angular  variable  is  used  among  the 
primary  DFs.  The  theoretical  (Hermann  and  Bagus  19771  and  experimental 
(Akimoto  et  al.  1979)  results  for  Ni(OOI)(ads)CO  support  this  assumption. 

The  remaining  problem  is  to  parametrize  the  off-diagonal  or  coupling  terms 
of  the  potential  energy  matrices,  labelled  by  an  index  "e".  Their  real  part 
represents  the  diabalic  electronic  coupling  between  the  basis  stales,  so  it  can 
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differ  significantly  from  zero  only  in  the  regions  close  to  crossings  or  avoided 
crossings  of  the  PESs.  These  regions  are  expected  to  occur  at  short  adsorbate- 
substrate  distances,  because  of  the  functional  form  of  the  Morse  potential.  One 
possible  form  to  describe  such  behaviour  is  with  a cut-off  function 

Vc(R,t)  = - tanh  [Qe(ft  - TT  - Re  + yrc)\)  (5.12) 

where  Re  and  rc  denote  the  position  of  surface  crossing. 

The  imaginary  coupling  potential  Wc(R,r)  appears,  as  was  deduced  in  Chapter 
3,  from  the  coupling  of  the  adsorbate  electronic  stales  with  the  interband  electronic 
transitions  (between  the  conduction  and  upper  empty  bands)  of  the  substrate.  Like 
Ws(R.r)  and  We(R,r),  it  is  expected  to  have  a short  range  contribution  coming  from 
the  choice  of  electronic  representation,  but  in  addition,  there  should  be  a longer 
range  term  from  the  interaction  of  the  induced  dipole  moments  of  the  adsorbate 
and  the  substrate.  The  proposed  parametrization  is 

H'o(R,r)  = ^y-{l  - tanli  (oc(/2  — *yr  — Rc  + 7Te))}  (5.13) 

which  comes  from  the  diagonalization  of  W and  is  analogous  to  the  real  diabatic 
coupling  Vi. 

The  real  PESs  (5.8)  and  (5.10)  with  diabatic  coupling  (5.12),  as  functions 
of  the  adsorbate-substrate  distance  R.  for  a fixed  value  of  r corresponding  to  the 
ground  electronic  state  equilibrium  vibrational  distance  between  C and  O.  arc 
illustrated  in  figure  5.1. 
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Figure  5.1:  Poicnli.il  energy  surfaces  with  diabolic  coupling  vs.  adsorbate-substrate 
distance  R for  fixed  equilibrium  value  t=r,.  for  the  Ni(OOI)(ads)CO  system. 


Transition  Dipole  Moment  for  CO  Photodesorption 


By  the  initial  assumption  of  this  study,  the  transition  dipole  moment  between 
the  adsorbed  electronic  state  of  CO  and  the  desorbing  excited  state  is  defined 
as  a function  of  primary  degrees  of  freedom  (i.e.  adsorbate  nuclear  coordinates) 

D = D(/l,r)  (5.14) 

The  customary  parametrizations  for  this  operator  in  the  literature  (Kirby  and 
Goldfield  1990.  Pradhan  et  a!.  1991)  arc  of  polynomial  type,  often  with  an 
attenuation  factor  to  account  for  the  relatively  short  range  of  interaction.  A more 
accurate  paramelrization  was  possible  due  to  the  electronic  structure  calculations, 
commented  in  the  next  section.  The  transition  dipole  was  found  to  increase  almost 
linearly  with  the  adsorbate-substrate  distance  R,  with  a gradual  cutoff  beyond  an 
R = R&  = const.-,  while  its  value  is  almost  independent  of  the  C-0  distance  r. 
again  up  to  a cutoff  at  r = r„(R).  A good  fitting  was  found  for  a function 
D(R,t)  = (D„  + D,n  + D,/i2){l  - tanli  |,(fi  - /^r)]}x 

(5.15) 

x{l-tanh(C(r-r„(«))]} 


Electronic  Structure  Calculations  of  the  Region  of  Primary  DFs 
This  section  gives  a more  detailed  explanation  of  the  electronic  structure 
calculations  that  were  done  for  the  NiCO  molecule  and  for  CO  adsorbed  on  the 
top  site  of  an  Nil0  cluster,  in  order  to  gain  a better  insight  into  the  CO  adsorption 
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mechanism  on  ihe  Ni(00l)  surface  and  ihe  character  of  the  electronic  transition 
that  leads  to  its  desorption,  and  to  obtain  a reliable  scmiempirical  parametrization 
for  the  PESs  with  their  couplings  and  the  transition  dipole  moment,  to  be  used 
later  in  the  compulation  of  the  primary  DF  dynamics.  These  calculations  were 
done  using  the  ZINDO  software  [Ridley  and  Zemer  1973,  Zemer  ct  al.  1980). 

To  calculate  the  ground  state  energy  of  the  given  molecule,  as  well  as  the 
effective  atomic  charges  and  electronic  populations  of  the  molecular  orbitals 
(MOs),  the  ZINDO  program  runs  an  iterative  Hartrec-Fock  (HF)  type  calculation 
for  a specified  geometry  and  total  spin  state.  The  relative  positions  of  the  C and  O 
atoms  with  respect  to  the  adsorption  site  Ni  atom  were  varied  around  the  known 
equilibrium  values  [Heskett  ct  al.  1985),  while  the  Ni10  cluster  was  constructed 
by  taking  into  account  the  eight  nearest  neighbours  and  the  atom  "underneath" 
the  adsorption  site  of  the  known  FCC  crystal  structure  of  the  nickel  solid.  The 
total  spin  state  was  chosen  to  be  a singlet  (no  unpaired  electrons)  for  NiCO  and 
a triplet  (two  unpaired  electrons)  for  Ni|0CO.  as  those  choices  were  the  closest 
approximations  to  the  experimental  value  of  the  effective  magnetic  moment  of  0.6 
Bohr  magnetons  per  atom  for  the  nickel  metal  (Kitlel  1957),  which  comes  from 
the  spin  decoupling  of  its  3d-4s  conducting  band.  In  fact,  a higher  multiplicity 
would  have  been  even  closer  for  NiioCO,  but  high  spin  multiplicities  present  a 
difficulty  for  convergence  of  the  calculations.  Information  on  the  input  data  for 
the  mentioned  ZINDO  calculations  can  be  found  in  Appendix  B. 

As  a general  trend,  the  HF  calculations  reproduced  quite  well  the  experimental 
equilibrium  O-C  and  C-Ni  distances  [Heskett  et  al.  1985);  however,  the  computed 
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binding  energy  was  substantially  Iatger  than  its  experimental  value,  due  to  the 
strong  deformation  of  the  MOs  of  CO  and  the  adsorption  site  Ni  atom.  This 
deformation,  in  fact  larger  for  Ni,0CO  'than  for  NiCO,  is  a consequence  of 
representing  an  infinite  solid  with  a small  atomic  cluster,  so  the  attractive  forces 
between  CO  and  the  Ni  atomic  cores  are  not  properly  balanced.  This  problem 
was  resolved  by  rescaling  the  computed  values  of  energy  to  fit  the  known  ground 
state  PES  (Andersson  and  Pendry  1979];  later,  energies  of  the  excited  states  were 
rescaled  by  the  same  factor  to  obtain  the  excited  PES  parametrization. 

As  staled  before,  the  parameters  of  the  dipole  interaction  and  the  coupling 
between  R and  rin  the  ground  state  PES  [eqn.  (5.8)]  were  determined  by  the 
ZINDO-HF  calculations;  they  arc  given  in  table  5.II. 


•Rlble  5.II:  Table  of  parameters  for  the  dipole  coupling  and  the  R-r  coupling  for  the 
Ni(OOIXads)CO  ground  electronic  state,  from  electronic  structure  calculations. 


Magnitude  in  atomic  units 


To  compute  the  excited  states  of  the  molecule  of  interest,  the  ZINDO  software 
uses  configuration  interaction  (Cl)  calculations,  initiated  by  a HF  ground  state,  but 
with  a changed  basis  set.  Again,  positions  of  nuclei  and  the  total  spin  state  must 
be  specified  for  each  run.  The  ZINDO-CI  calculations  give  as  a result  a multitude 
of  excitations  with  transition  energies,  transition  dipole  moments,  and  oscillator 
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strengths  that  are  related  to  the  probabilities  of  corresponding  transitions.  The 
input  data  for  those  calculations  are  listed  in  Appendix  B. 

The  ZINDO-CI  calculations  were  done  on  the  NiCO  molecule  to  deter- 
mine the  desorbing  excited  state  PES  and  the  corresponding  transition  dipole 
moment.  The  strongest  transition  was  detected  from  a linear  combination  of 
0.30|ATf  3d,, 0.951^134,)  and  0.95|IVi  3d*,.,,)  +0.30|AT  3 d„)  orbitals, 
to  the  0.88|rVi  4p.)  - 0.37|C  2i)  orbital  of  the  NiCO  complex.  The  excited  PES 
was  then  determined  by  adding  the  rescaled  transition  energy  to  the  ground  state 
value  for  different  interatomic  distances.  As  expected  from  the  theoretical  analy- 
sis, the  "destination"  orbital  is  antibonding  between  nickel  and  carbon,  so  it  does 
indeed  lead  to  desorption,  as  assumed  by  the  MGR  model.  The  parameters  of  the 
eqn.  (5.10)  were  then  fitted  to  reproduce  the  obtained  excited  PES  energies  as 
functions  of  R and  r\  their  values  are  given  in  table  5.III. 


Table  5.111:  Parameter  table  for  the  Ni(001)(ads)CO  excited  electronic  stale  PES. 


3dx>-,i)  - 0.95|Afi  3d,,)  and  0.95|(Vi  3dxi^)  + 0.30|/Vi  3d,,) 


lion  of0.30|Wi: 
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orbitals,  and  the  0.88|tVi  4p_-)  - 0.37|C  2a)  orbital  of  NiCO,  corresponding  to  the 
eqn  (5.15)  with  (5.16).  are  given  in  table  5.IV. 


CHAPTER  6 

NUMERICAL  ASPECTS  OF  TIME  PROPAGATION 


Methods  of  Time  Propagation  of  Wavcpackcts  with  Dissipation 
The  time  propagation  of  wavopackets  is  the  most  important  computational  as- 
pect of  this  study,  as  the  time  dependent  amplitudes  are  required  for  calculation  of 
the  time  correlation  functions  used  in  the  description  of  the  physics  of  photodes- 
orption.  This  propagation  in  fact  means  solving  numerically  the  lime  dependent 
Schrddinger  equation  (TDSE)  — several  methods  are  known  for  this  purpose  and 
can  be  found  in  literature  (Leforestier  ct  al.  1991,  Truong  et  al.  1991). 

The  problem  amounts  to  accurate  and  efficient  construction  of  the  time 
evolution  operator  of  the  TDSE,  which  for  the  primary  effective  Hamiltonian 

tf(A0=  exp  (-iAfP,,)  (6.1) 

with  h = 1.  This  operator  advances  the  solution  of  the  TDSE  by  the  given  lime 
step  At.  The  Hamiltonian  operator  Pp  has  been  derived  in  previous  chapters  and 

Fp  = Ppfi  ~ ti/dX ) = k^-iO/Ox')  + !/(*)  - iw(x)/l  (6,2) 

The  radiation  held  and  its  coupling  to  the  system  do  not  appear  in  this  Hamiltonian, 
which  implies  the  initial  conditions  of  the  time  evolution  such  that  the  light  had 
already  been  absorbed  by  the  system  and  the  primary  region  is  in  the  state  |e,0g) 
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(cf.  Chapter  4)  at  the  initial  time.  There  are  various  methods  of  numerical  time 
evolution  with  the  operator  (6.1)  [Leforesticr  el  al.  1991.  Truong  el  al.  1991); 
two  of  them  will  be  discussed  in  the  rest  of  this  section. 

The  first  one  starts  by  defining  a basis  set  of  Nb  suitable  functions 
[Hamilton  and  Light  1986]  that  spans  the  whole  space  of  interest  (in  this  context 
Nb  is  the  dimensionality  of  the  basis  set  and  not  the  physical  dimensionality  of 
space  — number  of  degrees  of  freedom).  The  wavepacket  0"(r)  then  becomes 
a column  matrix  with  elements 

*"{*)  = j i‘mr(X)dX  (6.3) 

and  the  kinetic  and  potential  energy  operator  become  matrices  with  elements 

Kjt  = ~Tm  (6.4) 

Vy»  = f 4){X)V{X)ik{X)iX 

r (6-5) 

Wyt  = J 4'(g)W(X)MX)dX 

so  that  the  Hamiltonian  matrix  is  the  sum  of  these  three  matrices 

H = K + V — iW/2  (6.6) 

Finally,  the  time  evolution  matrix,  allowing  for  a nonorthogonal  basis,  is 

U(A()  = exp  (— iAIS_1H)  (6.7) 

where  S is  the  overlap  matrix  of  the  basis  functions. 

The  advantage  of  the  basis  expansion  for  wavepackets  and  the  time  evolution 
operator  is  that  a good  choice  of  the  basis  set,  according  to  the  physical  properties 
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of  the  system  (symmetry,  etc.),  allows  the  use  of  relatively  few  basis  functions 
(Nj,  not  very  large)  and  thus  relatively  small  matrices  for  the  numerical  procedure. 
However,  Ng  does  become  prohibitively  large  for  relatively  extended  systems  with 
several  degrees  of  freedom.  Also,  the  matrix  elements  of  H are  integrals  that  must 
be  calculated;  and  H.  as  well  as  S for  a non-orthonormal  basis,  are  non-diagonal 
matrices,  so  they  must  be  diagonalized  to  construct  the  time  evolution  matrix. 

The  alternative  approach  is  a discretization  of  space  and  subsequent  evaluation 
of  the  Hamiltonian  matrix  elements  as  functional  values  at  the  given  points. 
The  space  of  primary  DF  coordinates  R and  r is  discretized  for  the  numerical 
implementation  of  the  time  evolution  scheme.  This  is  done  in  the  simplest 
fashion  by  introducing  an  equidistant  two-dimensional  discrete  grid  of  n (for 
the  R direction)  by  m (for  the  r direction)  points,  covering  the  whole  space  of 
interest;  the  discretized  variables  become  if,-  = Ro  + A R and  rj  = ro  + Ar, 
where  A R and  A r are  the  grid  spacings  in  R and  r,  respectively.  Although 
this  is  not  the  most  economical  way  in  terms  of  the  total  number  of  points 
necessary  (cf.  discrete  variable  representation  (Light  el  al.  1984.  Colbert  and 
Miller  1991,  Seidemann  and  Miller  1991)).  it  has  an  advantage  of  allowing  for 
different  numerical  algorithms,  such  as  fast  Fourier  transform  (FFT)  [Press  et  al. 
1992).  The  space  discretization  introduces  a coordinate  representation,  so  that 
(/ii|/i|/i;.)  = I'A/Ji,-  ;..  An  observable  — operator  A is  said  to  be  local  if 

= [A]id-d;y4jj.  (6.8) 

so  when  it  operates  on  a wavefunction  in  the  space  of  R and  r.  the  different 
components  in  the  coordinate  representation  do  not  mix.  Obviously,  all  operators 
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lhai  are  functions  of  coordinate  operators  P and  f are  local;  in  particular  the 
potential  energy  operators  in  both  diabatic  and  adiabatic  electronic  representation, 
as  well  as  in  the  "representation"  of  the  diagonal  imaginaty  potential.  Using  this 
property,  it  is  advantageous  to  use  the  split  operator  scheme,  defined  in  Chapter  4 
U<d>(A<)  = (M*  ■ (At/2)  - M]  - (L*  • U^(A'/2)  • LJ- 

,,,  , n , , («■’) 

•Ur  (At)  ■ [L1  • Uy  (At/2)  • L)  • (M*  ■ U$(Al/2)  ■ M) 

This  time  evolution  scheme  was  used  in  all  the  calculations  of  this  study. 

It  is  obvious  that  the  kinetic  energy  operator  k does  not  satisfy  the  locality 
condition  (6.8),  because  the  second  derivatives  with  respect  to  R and  r must 
be  expressed  c.g.  in  terms  of  finite  differences.  However,  k is  local  in  the 
momentum  representation,  as 


where  P and  fi  are  conjugate  momenta  to  R and  r.  respectively.  The  momentum 
space  is  discretized  according  to  Pk  = Pa  + kAP  and  « = po  + lAp.  where 
the  grid  spacings  arc  AP  = l/(nA fl)  and  Ap  = l/(mAr);  the  kinetic  energy 
operator  then  becomes 

k = El^.W)K[d)(Pl,Pt)(Pi,p,|  (6.1 1) 

where  K<d|  is  a diagonal  (nxm)xfnxml  matrix.  The  kinetic  energy  matrix 


j,  = E!/i'''Vl5,W)K'd>(a,w|/It.,ey.) 


(6.12) 


K<d>  = F^.  • K<d>  • Ft 


(613) 


of  the 


"(At)  = (l4  ■ U“( At/2)  • ]Lw  - (L^  ■ U<,*>(A</2)  - 1 
[F^U^AO-Ft]- 

l4 ' U(v*’(A‘/2)  • M • [L(v  ■ U(w‘(At/2)  - Lw] 
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Initial  Conditions  for  Wavcpackct  Propagation 
Initial  Conditions  for  the  Wavefunction 

As  stated  earlier,  the  primary  DF  space  of  CO  pholodesorption  from  the 
Ni(OOI)  surface  is  defined  in  terms  of  two  coordinates:  the  distance  between 
the  CO  center  of  mass  and  the  surface  R and  the  internal  C-0  distance  r.  This 
space  has  been  discretized  and  the  maximal  allowed  grid  spacings  are  generally 
determined  by  the  momentum  of  wavepackcts  that  propagate  over  it  in  each 
direction.  The  initial  state  of  the  system  is  described  by  a two-component 
wavefunction 


This  initial  slate  means  that  the  time  evolution  begins  after  a short  light  pulse  has 
been  shining  on  the  adsorbate-substrate  complex,  photons  have  been  absorbed  and 
a transition  into  the  excited  electronic  state  of  the  adsorbate-substrate  complex 
has  occurred. 

The  exact  form  of  the  wavcpacket  0)  depends  on  the  exact  charac- 

teristics of  the  electronic  transition.  In  general,  it  is  assumed  that  the  transition 
has  a Franck-Condon  character,  i.e.  that  the  expectation  value  of  the  two  coor- 
dinates and  their  conjugate  momenta,  as  well  as  the  shape  of  the  wavcpacket, 
remain  unchanged.  Therefore,  if  the  system  had  originally  been  in  its  ground 
vibrational  state  in  the  Morse-harmonic  potential  well  of  the  ground  electronic 
PES,  tfe(/2,r;  0)  can  be  approximated  by  a two-dimensional  harmonic  oscillator 
ground  stale  eigenfunction  — a Gaussian  in  both  R and  r,  i.e. 

MR. no)  = Wexp  [-4(72  - flo)*]  exp  [-o2(r  - r0)2] 


(6.17) 
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where  S0  and  ro  denote  the  initial  position,  the  width  parameters  are  a,  = 
Vscofis/A  and  qr  is  obtained  from  the  shape  of  the  Morse  potential  well  of  the 
ground  stale  PES,  and  N is  a normalization  factor.  This  is  certainly  trac  in  the 
case  of  resonant  absorption,  when  the  energy  of  the  absorbed  photon  is  exactly 
equal  to  the  energy  difference  between  the  two  PESs  at  the  initial  position  of  the 
adsorbate,  reduced  by  the  zero-level  vibrational  energy  in  /?. 

If  the  energy  of  the  absorbed  photon  is  above  the  resonant  absorption  value, 
the  initial  wavepacket  must  have  an  average  kinetic  energy  equal  to  the  excess 
of  the  photon  energy,  even  though  its  average  momentum  must  still  vanish.  One 
way  to  construct  such  a wavepacket  is  to  multiply  the  original  packet  by  a phase 
factor  exp  ( iRP ) and  then  integrate  over  a distribution  of  momenta  with  a standard 
deviation  such  that  (P2)  - 2Afcoh  Aw,  where  Aw  is  the  difference  of  the  photon 
frequency  from  the  resonant  absorption  value.  The  resulting  wavepacket  has  a 


= N exp  [-(aj  + MCoAw/A)(/!  - flo)2]  exp  [-ar2(r  - r0)2] 
(6.18) 


in  the  coordinate  space  — it  becomes  narrower  then  the  original  wavepacket  in  R. 
This  should  be  acceptable  for  the  values  of  Aw  small  compared  to  the  vibrational 
quanta  in  the  excited  electronic  state,  but  it  can  lead  to  unphysical  results  when 
Aw  becomes  relatively  large.  In  that  case  it  may  be  expected  that  the  CO  internal 
vibrational  levels  will  become  populated,  instead  of  converting  the  whole  excess 
of  the  photon  energy  into  the  kineUc  energy  of  the  desorbing  species. 


Initial  Conditions  for  the  Dipole  Amplitude 


The  initial  conditions  described  above  refer  to  the  propagation  of  the  wave* 
function,  which  is  suitable  to  study  the  experimental  situations  involving  excita- 
tions  by  short  light  pulses.  If.  on  the  other  hand,  one  wishes  to  study  a steady  state 
experimental  situation  in  which  the  adsorbate-substrate  complex  interacts  with  a 
steady  beam  of  light,  the  physically  appropriate  description  is  the  propagation  of 
the  dipole  amplitude  and  the  above  initial  wavcpackct  must  be  multiplied  by  the 
transition  dipole  moment.  In  that  case,  the  time  evolution  operator  becomes  an 
effective  frequency  operator  (cf.  Chapter  2)  [Srivastava  and  Micha  1 1991)  — the 
initial  energy  of  the  system  is  subtracted  from  the  primary  effective  Hamiltonian. 

= exp  [-i(Fp  - eJ  t/hj  (6.19) 

The  initial  wavcpacket  in  this  case  must  represent  the  initial  transition  dipole 
amplitude,  defined  in  Chapter  2;  therefore. 

and  r;0)  is  the  initial  wavcfunction  defined  by  the  eqn.  (6.17).  In  this 
case,  there  arc  no  modifications  coming  from  Ihc  frequency  of  the  absorbed  light; 
in  fact,  the  Fourier  transform  with  respect  to  time  of  the  dipole-dipole  TCF  is  a 
function  of  the  photon  frequency  and  it  gives  the  probability  of  desorption  after 
absorption  of  a photon  of  the  given  frequency. 


Numerical  Parameters  for  Time  Propagaiic 


Al  later  limes  r > 0,  the  wavcpacket  will  have  an  average  momentum  of 
(P)t  = £'  E (MOI-MO)  »nd  similar  for  the  r 

direction.  The  maximal  allowable  grid  spacings  are  then  estimated  by 


AR  « 2t/(P), 


(6.21) 


At  « 2x /(p), 

where  "much  less  than”  implies  at  least  several  grid  points  per  unit  wavelength, 
i.e.  about  an  order  of  magnitude. 

The  other  critical  parameter  of  the  numerical  time  propagation  is  the  time  step. 
Neglecting  for  the  moment  the  potentials,  the  time  evolution  operator  matrix  on 
a grid  with  spacing  AR  for  a free  particle  of  mass  M can  be  written  as 


U(At)=exp  -i-  .D 

[ 2M(ARf  \ 

where  D is  a matrix  of  finite  differences,  i.e.  its  elements  arc 
order  of  unity  and  smaller.  For  stability  of  propagation,  the  ra 
above  must  be  kept  relatively  small,  as  Uf Al)  is  multiplied  intt 
varies  in  space  according  to  the  momentum  of  the  particle. 


(6.22) 

all  numbers  of  the 
tio  in  the  exponent 
> an  amplitude  that 


The  time  evolution  method  developed  in  this  study  can  examine  not  only  the 
desorption  from  surfaces,  but  molecular  scattering  problems  as  well.  Test  calcu- 
lations were  performed  for  wavcpackcls  of  initial  Gaussian  amplitude,  modulated 
by  a phase  factor  exp  (-ifijfi)  to  reproduce  the  initial  momentum  in  the  molec- 
ular scattering  calculations.  The  method  was  tested  to  reproduce  the  results  of 
two  model  one-dimensional  (fixed  r)  scattering  problems  with  two  coupled  PESs 
and  without  dissipation  [Alvarcllos  and  Mctiu  1988.  Cohen  and  Micha  1992). 
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The  results  for  average  positions  and  electronic  slate  populations  as  functions  of 
time  [Alvarellos  and  Metiu  1988]  were  well  reproduced;  while  the  results  for  the 
electronic  transition  probability  as  a function  of  energy  were  in  agreement  with 
the  reference  [Cohen  and  Micha  1992)  at  low  energies  (less  than  0.4  Hartree),  the 
high  projectile  momenta  presented  a problem  in  the  time-dependent  calculations 
for  comparisons  at  the  higher  energies. 


The  typical  momenta  for  the  photodcsorption  problem  are  smaller  then  the 
upper  limit  that  was  tested  (for  example,  a CO  molecule  with  momentum  of  30  au 
has  the  kinetic  energy  of  approx.  0.25  eV)  so  the  calculations  could  be  performed 
with  larger  time  steps,  as  well  as  with  fewer  grid  points  in  each  direction. 


table  6.1  gives  the  numerical  propagation  parameters  in  both  the  preliminary  test 
calculations  and  the  actual  wavcpacket  propagation  used  in  calculations  of  the 
CO  photodesorption  from  nickel. 

Algorithm  for  the  Numerical  Time  Propagation 

The  numerical  propagation  of  wavepackcts  used  in  the  calculations  of  the  CO 
photodesorption  from  the  Ni(OOI)  surface  was  executed  with  a Fortran  77  program 
“cphe2d.r  and  the  corresponding  subroutines;  the  flow  diagram  of  this  algorithm 
is  given  in  Appendix  C.  The  algorithm  propagates  a two  dimensional  wavcpacket, 
given  as  a two  index  array  by  the  eqn.  (6. 15),  with  a time  evolution  operator  (6. 14) 
constructed  from  a Hamiltonian  operator  with  two  coupled  complex  potential 
energy  surfaces  by  means  of  a split  operator  scheme,  modified  to  include  the 
dissipation. 

The  algorithm  begins  by  the  input  of  data,  separated  in  two  subroutines  that 
read  in  the  data  from  two  separate  input  files.  The  first  one.  called  "inwh2d.r', 
reads  in  the  data  for  the  initial  wavefunction  given  by  the  eqn.  (6.17)  and  the 
transition  dipole  moment  with  a functional  form  discussed  in  Chapter  5.  The 
second  one,  named  “icpi2d.f\  reads  in  the  data  for  the  complex  PESs  and  their 
couplings,  all  defined  in  Chapter  5.  The  program  then  proceeds  by  constructing 
the  initial  wavcpacket  (subroutine  "inho2d.r).  and.  if  the  dipole-dipole  TCF 
is  to  be  calculated,  the  initial  dipole  amplitude  (subroutine  "dpam2d.f');  and 
also  the  complex  PESs  with  couplings  (subroutine  "cpts2dj”).  Evidently,  this 
structure  assures  a great  flexibility  in  the  use  of  the  program,  as  different  initial 
conditions  or  potential  energy  functions  can  be  obtained  by  merely  modifying 
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•he  corresponding  subroutines,  without  changing  the  main  program.  Before 
starting  the  time  propagation,  the  program  computes  the  initial  values  of  the 
time  dependent  properties  such  as  the  expectation  values  of  coordinate  operators 
A and  r (subroutine  "cent2d.f"),  and  the  TCFs  — autocorrelation  function  or 
dipole-dipole  TCF  (subroutine  “acrf2d.f"). 

The  time  evolution  in  this  algorithm  is  sequential  in  time  steps,  using  the 
modified  split  operator  scheme  (6.14)  to  advance  by  one  time  step  and  calculating 
the  time  dependent  properties  at  each  time  step  until  reaching  the  given  final 
time.  The  time  advancement  is  executed  as  a sequence  of  matrix  multiplications 
given  by  (6.14),  with  overwriting  of  the  wavepacket  after  each  multiplication;  for 
instance,  operating  with  the  imaginary  potential  propagator  is  done  as  a sequence 
X W -Lw-#,d)(l) 

x(l,  At)  =:Uwd(A(/2)- !*.<“>(!) 

...  (6-23) 

tA  d (l,  At)  =:  x(t,  At) 

X(t,At)  =:L(vi/i<d'(t,At) 

V><d>(t.At)  =:X(t,At) 

using  subroutines  "repi2d.f'  for  multiplications  with  L\y  and  its  adjoint,  and 
“adcp2d.r  for  multiplication  with  the  propagator  UwofAr/2).  An  analogous  pro- 


Mrepr2d.r'  and  "adpo2d.r.  The  time  advancement  by  the  kinetic  energy  propaga- 
tor is  done  with  a subroutine  "adft2d.r\  which  arranges  the  input  wavepacket  into 
a one  dimensional  array  suitable  far  the  fast  Fourier  transform  (FFT)  algorithm, 
takes  an  FFT  into  the  momentum  space  using  the  well  known  "fourn.r  routine 


(Press  ei  ai.  1992),  advances  by  a slep  A I in  the  momentum  space,  lakes  the 
inverse  FFT  back  into  the  coordinate  space  and  rearranges  the  output  wavepacket 
into  a two  dimensional  array,  overwriting  the  input  array. 

The  output  of  the  time  evolution  program  consists  of  the  time  dependent 
properties  — expectation  values  of  the  primary  coordinates.  TCFs  and  relative 
populations  of  the  two  electronic  states  (normalized  to  unity  at  each  time  step; 
the  actual  populations,  obtained  as  absolute  value  squares  of  each  wavefunction 
component  at  each  time,  go  asymptotically  to  very  small  limits  due  to  the 
dissipation),  as  well  as  the  input  files  to  obtain  the  two  dimensional  plots  of 
the  PESs  and  their  couplings,  the  transition  dipole  moment  between  the  electronic 
basis  states,  and  the  propagating  amplitudes  as  they  evolve  over  time. 

The  time  propagation  algorithm  described  in  this  section  achieves  good 
efficiency  in  both  computational  time  and  memory  occupation.  Although  the 
complex  PESs  and  their  couplings  are  stored  as  two  dimensional  arrays,  they 
need  not  be  generated  at  each  time  step  of  the  propagation  — this,  though 
increasing  the  necessary  memory,  considerably  accelerates  the  computation.  To 
illustrate  the  performance,  it  is  worth  noting  that  one  photodesorption  calculation 
with  parameters  from  the  third  column  of  the  table  6.1,  executed  on  an  IBM 
RS/6000-580  server  machine,  takes  about  106  minutes  of  CPU  lime,  occupying 


approximaiely  1.7  Mb  of  memory. 


CHAPTER  7 

RESULTS  OF  CALCULATIONS  AND  DISCUSSION 


Model  Calculations  for  Dissipative  Potential  Parameters 


This  section  addresses  the  parametrization  of  the  dissipative  imaginary  poten- 
tials that  come  from  the  SCF  treatment  of  the  adsorbate-substrate  coupling.  This 
problem  can  be  related  to  the  calculated  finite  lifetimes  of  the  CO  vibrational 
states  [Persson  and  Pcrsson  1980,  Head-Cordon  and  Tully  1992]  when  adsorbed 


on  metal  surfaces.  The  results  suggest  that  the  lifetimes  of  the  C-0  stretch  vibra- 
tions must  be  of  the  order  of  1 ps.  while  the  lifetimes  of  the  frustrated  translation 
mode  (C-Ni  stretch)  are  of  the  order  of  10  to  100  ps.  To  model  this  system,  a 
simplified  ground  state  complex  PES  is  defined  by 

W')  = - R,)1  + i-4(r  - rff)2— 

W.  (7I) 

where  the  frustrated  translation  frequency  fir,  is  obtained  by  expanding  the  ground 
state  Morse  potential,  given  by  the  cqn.  (5.1),  around  its  minimum.  No  coupling 
of  motions  in  R and  r was  included  because  it  was  found  to  be  small  (cf.  Chapter 
3)- 


The  harmonic  oscillator  eigenfunctions  in  two  dimensions  with  different 
quantum  numbers  (combinations  of  0,  I.  and  2)  were  then  propagated  with  the 
potential  (7.1)  in  order  to  test  the  computer  program  for  propagation  and  to  obtain 
an  estimate  on  the  dissipative  potential  strength  in  relation  to  the  lifetime.  The 


Chapter  2,  was  calculated  for  each 


87 

ihe  absence  of  dissipaiion,  its  absolute  value  square  is  equal  to  one  at  all  times, 
which  confirmed  that  the  propagation  algorithm  contained  no  spurious  couplings 
between  the  two  primary  coordinates.  An  example  — the  ground  vibrational  state 


in  figure  7.1. 

Assuming  an  exponential  decay  of  a given  vibrational  state,  the  lifetime  can 
be  expressed  as 


T"Wi'2,{~Sln[|c(1)li]}  (7-2) 

where  in  practice  the  limit  is  evaluated  for  the  final  time  of  the  propagation. 
'(/=!. Ox  I03  au.  Results  obtained  for  the  ground  harmonic  state  (0.0)  are  presented 
in  table  7.1 


Table  7.1:  Lifetime  of  ground  vibrational  stale  of  Ni(ads)CO 
for  different  values  of  the  imaginary  potential  strength  Wo. 


Wo  (10-3  Eh) 

1 C (If)  lJ 

r (ps) 

0.00 

1.0 

CO 

0.02 

0.97980 

1.0717 

0.05 

0.95031 

0.4291 

0.10 

0.90317 

0.2147 

0.20 

0.81598 

0.1071 

0.50 

0.60281 

0.0430 

1.00 

0.36609 

0.0218 
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Figure  7.1:  Autocorrelation  Function  vs.  time  for  a 2D  harmonic  oscillator 
ground  state,  for  different  values  of  the  dissipation  strength  parameter. 
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These  results  clearly  show  that  for  the  given  model  the  lifetime  of  vibrational 
states  ts  inverse  proportional  to  the  strength  of  the  dissipative  imaginary  potential; 
similar  results  were  obtained  for  the  excited  vibrational  states  of  the  system.  They 
indicate  that  the  imaginary  part  of  the  primary  ground  electronic  state  PES  must 
have  a magnitude  of  the  order  of  10-5  au  in  order  to  reproduce  the  known 
vibrational  lifetimes.  For  the  case  of  Ni(OOI)(ads)CO,  however,  the  electronic 
excitations  are  expected  to  be  the  primary  way  of  energy  dissipation  (Newns, 
Heinz  and  Miscwich  1991).  so  the  parameters  Vfi,  and  W°,  are  expected  to 
be  larger;  the  initial  guess  for  the  time  propagation  in  this  work  was  W°  = 
8.78  x 10  4 au  and  W®  = 1.14  x 10~3  au.  The  magnitude  of  the  imaginary  part 
of  the  diabatic  coupling  was  initially  set  at  W?  = 1.23  x 10-3  au. 

The  calculations  that  follow  were  done  for  two  different  initial  conditions, 
corresponding  to  an  excitation  by  a light  pulse  or  by  a steady  light  beam,  as 
indicated  in  Chapter  6. 

Calculations  of  CO  Photodesorption:  Varying  the  Incident  Light  Wavelength 

The  theory  of  photodesorption  dynamics  developed  in  this  study,  with  the 
initial  conditions  specified  in  Chapter  6,  assumes  up  to  this  point  that  the  electronic 
transition  in  the  adsorbate-substrate  complex  leading  to  the  desorption  was  caused 
by  a resonant  absorption  of  a photon,  i.e.  that  its  frequency  was  exactly  equal 
to  the  difference  of  potential  energies  of  the  two  electronic  states.  For  the 
PESs  parametrized  as  in  Chapter  5,  this  resonant  absorption  energy  is  = 
1.90  x 10~‘  au,  or  5.16  eV.  If  the  system  absorbs  a photon  with  a frequency  above 
the  resonance  value,  the  excess  of  energy  is  expected  to  be  partly  convened  into 


Ihe  kinetic  energy  of  rhe  adsorbate  as  it  moves  on  the  excited  PES,  and  partly 
dissipated  into  the  substrate.  This  section  presents  and  discusses  the  results  of  the 
dynamics  of  the  desorbing  species,  depending  upon  the  absorbed  photon  frequency 


Figure  7.2  shows  the  expectation  value  of  the  adsorbate-substrate  distance 
(«)r  = W'W|fl|^(0)/W)I^M>  » a function  of  time  for  four  different 
frequencies  of  the  absorbed  light.  The  initial  conditions  were  as  discussed  in 
Chapter  6,  for  excitation  by  a light  pulse.  As  expected,  the  asymptotic  velocity  of 
the  adsorbate  increases  with  the  frequency  shift  Aw;  it  is  possible  to  estimate 
it  directly  from  the  graph,  and  therefrom  calculate  the  kinetic  energy  of  the 
adsorbate  relative  to  the  substrate.  Comparing  the  shift  in  the  kinetic  cneigy 
AEk  = - Ek  (wj“)  with  the  frequency  shift  gives  information  about  the 

fraction  of  the  incident  photon  energy  that  is  transferred  to  the  relative  motion  of 
the  desorbing  species.  The  results,  given  in  the  table  7.II,  show  that,  although  A EK 
increases  monotonously  with  Aw,  the  increments  are  not  equal  nor  do  they  change 
at  the  same  rate.  The  reason  is  that  a change  in  the  absorption  frequency  and  a 
consequent  change  in  the  kinetic  energy  deform  the  primary  nuclear  wavefunction, 
and  the  energy  dissipation  in  the  substrate  is  altered  as  well.  This  is  apparent 
from  the  short  time  behaviour  of  the  average  ft  — as  Aw  increases,  it  appears  to 
decrease  for  short  time  as  some  components  of  the  wavepacket  acquire  negative 
values  of  the  relative  momentum;  later,  they  "rebound"  off  the  potential  barrier 
of  the  desorbing  state  and  the  wavefunction  moves  toward  large  ft  faster. 


Table  7.11:  Change  in  the  CO  kinetic  energy  as  a function  of  the  incident  photon  energy. 


Figure  7.3  shows  the  absolute  value  square  of  the  autocorrelation  function, 
followed  by  the  plot  of  its  natural  logarithm,  given  by  figure  7.4.  It  is  apparent  that 
the  autocorrelation  function  decays  more  slowly  at  short  times  if  the  absorption 
frequency  is  increased,  which  can  again  be  attributed  to  the  deformation  of  the 
wavefunction  due  to  the  increasing  kinetic  energy.  Figure  7.4  shows  that  the 
late  time  bohaviour  of  the  autocorrelation  function  does  not  change  significantly 
within  the  range  through  which  the  absorption  frequency  was  varied  in  the 
calculations.  The  desorption  lifetime  can  then  be  deduced  from  the  cqn.  (7.2),  it 
is  approximately  equal  to  0.16  picoseconds  for  the  case  of  resonant  absorption. 

Having  used  initial  conditions  for  excitation  by  a steady  light  beam,  figure  7.3 
shows  the  absolute  value  square  of  the  dipole-dipole  TCF  for  the  Ni(OOI)(ads)CO 
system,  for  the  varying  frequency  of  incident  light;  the  transition  dipole  moment 
parametrization  from  Chapter  5 was  used  in  the  construction  of  the  initial  dipole 
amplitude.  It  is  evident  that  the  increase  in  the  absorbed  photon  frequency  causes 
broadening  and  a stronger  pronunciation  of  the  structure  in  the  dipole-dipole  TCF. 
The  desorption  lineshape  is  obtained,  as  shown  in  Chapter  2,  as  a Fourier  transform 


are  shown  in  figure  7.6. 
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Figure  7.2:  Average  adsorbate-substrate  distance  R vs.  lime  for  different  values  of  the 
absorbed  photon  frequency;  Aui  is  the  shift  from  the  resonant  absorption  frequency. 


Figure  7.3:  Absolute  value  square  of  the  autocorrelation  function 
vs.  time  for  different  values  of  the  absorbed  photon  frequency. 
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lime  (au) 


Figure  7.5:  Absolute  value  square  of  the  dipole-dipole  time  correlation 
function  vs,  time  for  different  values  of  the  absorbed  photon  frequency. 


CO(eV) 


Figure  7.6;  Pholodesorption  lmcshape  vs.  frequency 
for  different  values  of  the  absorbed  photon  frequency. 
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The  iincshapc  shows  a peak  at  a frequency  close  to  the  frequency  of  the 
absorbed  light,  and  the  position  of  the  peak  changes  in  agreement  with  the 
absorption  frequency  increase,  though  always  remaining  somewhat  lower.  There 
are  two  reasons  for  this:  the  presence  of  the  dissipation  and  the  paramctrization  of 
the  PESs  itself — the  excitation  energy  was  not  adjusted  to  the  experimental  value 
and  the  calculated  value  was  retained  instead.  There  is  a well  detectable  secondary 
maximum  in  the  lineshape,  approximately  1 eV  above  the  peak  frequency;  it  is 
attributable  to  the  internal  vibrational  excitations  in  the  adsorbate. 

Calculations  of  CO  Photodcsorption:  Varying  the  Dissipation  Strength 

Developing  the  physical  model  of  photodcsorption,  it  was  shown  in  Chapter 
3 that  the  magnitude  of  the  dissipative  imaginary  potentials  is  directly  related 
to  both  the  pair  correiauon  functions  of  the  lattice  displacements  and  of  the 
electronic  charge  density  in  the  substrate.  However,  the  detailed  dynamics  of 
the  metal  substrate  is  beyond  the  scope  of  this  study,  so  the  dissipative  potentials 
were  parametrized  semiempirically  as  the  imaginary  parts  of  the  diabatic  PESs 
and  their  coupling  (cf.  Chapter  5).  They  were  chosen  to  be  decaying  exponential 
functions  (for  diagonal  terms  in  the  electronic  basis),  and  a switching  function 
(for  the  diabatic  coupling  term);  with  the  same  inverse  range  parameters  as  the 
corresponding  potential  energy  functions.  The  magnitudes  of  the  dissipation 
potentials  and  coupling  could  not  be  determined  by  the  electronic  structure 
calculations;  the  model  calculations  from  the  beginning  of  this  Chapter  have 
provided  only  a qualitative  insight  to  those  magnitudes,  by  comparison  with  the 
known  lifetimes  of  vibrationally  excited  states  of  the  adsorbate.  This  section 


presents  the  results  of  the  CO  photodesorption  dynamics  as  they  depend  on  the 
magnitude  of  the  dissipative  potentials  and  their  coupling.  The  corresponding 
parameters  W,°.  W,°  and  Wt°  were  varied  from  the  values  coming  from  the 
model  calculations  to  values  2.  4 and  10  times  larger. 


First,  the  time  dependence  of  the  expectation  values  for  the  primary  coordi- 
nates R and  r is  shown  for  different  values  of  the  dissipation  parameters  on  figures 
7.7  and  7.8.  respectively.  The  qualitative  behaviour  of  both  coordinates  shows  up 
as  expected  — increasing  the  dissipation  strength  lowers  the  asymptotic  velocity 
of  the  desorbing  species  and  also  reduces  the  amplitude  of  its  internal  vibration. 


To  examine  more  precisely  the  dependence  of  the  dynamics  upon 
strength,  one  can  calculate,  according  to  the  classical  corresponder 
energy  of  the  desorbing  species 


Is  internal  vibrational  energy 
r _ KoSc  r,_\2 


the  dissipation 
ice.  the  kinetic 

03) 

(7.4) 


where  (r)m„  denotes  the  maximum  displacement  in  r from  its  equilibrium 
position.  It  must  be  noted  here  that  the  classical  result  for  the  vibrational  energy 
cannot  be  expected  to  give  an  accurate  result  as  it  does  not  account  for  the  zero 
point  energy;  still,  it  can  illustrate  the  effect  of  dissipation  qualitatively.  The 
results  for  £»■  and  £,  arc  summarized  in  table  7.III,  they  show  that  the  reduction 
in  the  kinetic  and  the  internal  vibrational  energy  of  the  adsorbate  is  approximately 
proportional  to  the  increase  in  the  dissipation  strength. 
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Figure  7.7:  Expectation  value  of  the  adsorbate-substrate 
distance  R vs.  time  for  varying  dissipation  strength. 
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Analysis  of  the  energy  dissipation  effects  on  the  photodesorption  process 
continues  by  estimating  the  lifetime  of  desorbing  states  for  different  values  of 
the  dissipation  strength  parameters.  As  before,  the  lifetime  is  deduced  from  the 
asymptotic  slope  of  the  logarithm  of  the  absolute  value  square  of  the  autocorrela- 
tion function,  which  is  shown  in  figure  7.9  for  the  dissipation  strength  parameters 
indicated  before.  The  results  for  the  photodesorption  lifetime  are  given  in  table 
7JV;  they  clearly  show  that  the  desorption  lifetime  is  inverse  proportional  to  the 
dissipation  strength,  as  in  the  model  calculations  at  the  beginning  of  the  Chapter. 


Dissipation  strength  parr 


It  is  also  of  interest  to  examine  the  effects  of  the  energy  dissipation  in  the 
substrate  on  a desorbing  species  that  had  initially  been  in  a stationary  state,  after  an 
electronic  transition  caused  by  a stationary  light  beam  at  the  resonant  absorption 
frequency.  To  this  effect,  the  dipole-dipole  TCF  (figure  7.10)  and  the  desotplion 
lineshape  (figure  7.1 1)  were  calculated  for  each  set  of  the  dissipation  strength 


parameters. 


time  (au) 
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Figure  7.10:  Absolute  value  square  of  Ihe  dipole-dipole  TCF  vs,  lime  for 
different  values  of  the  dissipation  strength,  after  a resonant  absorption  of  light. 


Figure  7.11:  Phoiodcsorpuon  lincshape  vs.  frequency  for  different 
values  of  the  dissipation  sirengih.  after  a resonam  absorption  of  lighi. 
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As  could  be  expected,  the  increase  in  the  dissipation  strength  causes  a faster 
attenu.»ion  of  the  dipole-dipole  TCF  over  time.  As  a consequence,  the  lineshapc's 
largest  peak,  obtained  from  the  Fourier  transform  of  the  TCF,  broadens.  It  is  also 
notable  that  the  relative  intensity  of  the  secondary  maximum  increases,  which 
means  that  the  stronger  dissipative  coupling  to  the  substrate  causes  a stronger 
internal  vibrational  excitation  in  the  adsorbate. 

The  dissipative  coupling  between  the  electronic  states  of  the  adsorbate- 
substrate  complex  accounts  for  the  possibility  of  electronic  transitions  between 
the  two  electronic  basis  states  — those  transitions  are  mainly  the  deexcitation 
of  the  Ni  substrate  from  the  4p  empty  band  slates  to  the  3d/4s  conduction  band 
states.  Their  occurrence  is  manifested  as  the  relative  population  of  the  ground 
electronic  stale,  defined  as 

m = W‘>K(<)>/KtfJ«W0>  + «(<M(0»  (7.5) 

This  relative  population  is  shown  in  figure  7. 12  for  the  varying  dissipation  strength 
parameters.  The  results  show  that  in  all  cases  an  electronic  transition  to  the  ground 
stale  PES  is  more  probable  and  will  occur  at  shorter  times;  when  the  dissipative 
coupling  is  increased,  the  short  time  deexcitation  is  more  pronounced,  but  it 
disappears  sooner  at  long  times.  It  is  apparent  from  figure  7.12  that  the  maximal 
probability  of  the  electronic  deexcitation  increases  linearly  with  the  dissipative 
coupling  strength. 


Figure  7.12:  Relative  population  on  the  ground  state 
PES  vs.  time,  for  varying  dissipation  strength. 


Calculations  of  CO  Photodcsorption:  Varying  the 
Asymptotic  Energy  Difference  between  PESs 

In  this  study’s  model  of  the  CO  photodesorption  from  nickel,  the  product 
states  consist  of  the  ground  electronic  state  of  the  desorbing  CO  with  an  internal 
vibrational  distribution,  and  a 4p  or  a 3d/4s  Ni  electronic  band  state.  The  energy 
difference  between  the  corresponding  PESs  at  large  adsorbate-substrate  distances 
then  represents  the  energy  difference  between  the  two  metal  bands;  the  value  of 
AE  = 3.48  eV  was  obtained  by  extrapolating  the  results  of  the  NiCO  molecule 
electronic  structure  calculations  to  large  R with  an  exponentially  decaying  function 
for  the  excited  state  PES,  and  it  coincides  with  the  energy  difference  between  the 
Fermi  level  and  the  maximum  of  the  4p  band  density  of  states  [Nagel  1973J  for 
solid  nickel.  However,  since  electronic  band  slates  form  an  energy  continuum,  the 
asymptotic  energy  value  on  the  excited  state  PES  cannot  be  fixed  with  certainty. 
Therefore,  the  parameter  AE  of  the  potential  energy  function  V,(R.rj  (cf.  Chapter 
5)  was  changed  in  the  calculations  by  adding  a long  range  correction  of  the  form 


v“”(H,r)  = Vc(R,t)  - l£,(l  + truth  [Q«(fl  - 2/iaI))}/2  (7.6) 

which  lowers  the  asymptotic  energy  difference.  The  calculations  were  repeated 
for  AE  = 2.98  eV,  and  both  the  autocorrelation  and  the  dipole-dipole  TCF 
were  virtually  unchanged;  thus  the  0.5  eV  decrease  in  the  asymptotic  energy 
difference  was  not  found  to  affect  with  any  significance  the  photodesorption 
lifetime  and  lineshape.  The  only  change  was  shown  by  the  late  time  behaviour 
of  the  expectation  value  of  R,  illustrated  by  figure  7.13;  however,  the  calculated 


109 

increase  in  ihc  final  kineiic  energy  of  ihc  desorbing  species  is  only  AEk  = 0.054 
cV  or  1 1%  of  ihc  change  in  A E. 
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Figure  7.13:  Average  adsorbate-substrate  distance  R vs. 
time,  varying  the  asymptotic  energy  difference  AE. 


CHAPTER  8 
CONCLUSIONS 

Contributions  to  the  Theory  of  Photodcsorption 

The  physical  model  developed  in  this  dissertation  describes  the  desorption 
induced  by  an  electronic  transition  (DIET),  caused  by  an  interaction  of  UV  light 
with  a localized  region  of  the  adsorbate-substrate  complex.  This  establishes  a 
hierarchy  of  the  degrees  of  freedom  (DFs)  in  the  complex:  those  that  involve 
the  operator  of  the  transition  into  a desorbing  state  are  the  primary  DFs  and 
they  describe  the  relative  motion  of  the  desorbing  species  with  respect  to  the 
substrate  as  well  as  its  internal  motion,  others  arc  the  secondary  DFs  and  they 
describe  the  substrate  dynamics.  Primary  and  secondary  DFs  are  coupled  sclf- 
consistently.  which  leads  to  the  energy  fluctuation  and  dissipation  in  the  primary 
DF  Hamiltonian  of  the  desorbing  species.  Assuming  a stochastic  character  for 
the  secondary  DF  dynamics,  the  energy  fluctuations  vanish  on  the  average  and 
the  dissipation  occurs  through  relaxation  into  a temporary  equilibrium  of  the 
secondary  DFs,  corresponding  to  a time  dependent  temperature  of  the  substrate. 

The  energy  dissipation  term  appears  in  the  primary  DF  Hamiltonian  as  an  anti- 
Hcrmitian  addition  to  the  potential  energy  function,  resembling  the  imaginary 
optical  potentials  from  scattering  theory.  Its  form  was  derived  for  two  types 
of  substrate  excitations:  lattice  vibrations  — phonons  and  electron-hole  pair 
excitations.  In  both  cases,  the  form  of  the  dissipative  potential  depends  upon 
a time  correlation  function  (TCF)  of  a substrate  property.  For  phonons,  the 
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magnitude  of  dissipation  is  proportional  to  the  TCF  of  the  lattice  displacements; 
while  for  the  electronic  excitations,  it  is  proportional  to  the  charge  density 
fluctuation  TCF.  The  self-consistent  field  (SCF)  coupling  of  the  primary  and 
secondary  DFs  is  crucial  to  the  dissipative  potential  as  both  the  displacement  and 
the  charge  density  fluctuation  TCFs  vanish  in  its  absence;  therefore,  the  energy 
d.M.Ka..v..  m the  desorbing  species  dynamics  is  a direct  consequence  of  the  SCF 
treatment  of  the  system  that  was  developed  in  this  work. 

The  experimental  conditions  of  pholodesorption  with  electronic  transitions 
may  involve  either  the  adsorbate  interacting  with  a short  pulse  of  light  inducing 
an  approximately  instantaneous  transition  to  the  desorbing  state,  or  a steady  beam 
of  light  preparing  a stationary  initial  state  for  desorption.  In  the  former  case,  it 
is  appropriate  to  measure  the  lifetime  of  the  desorbing  quantum  state  after  the 
transition;  in  the  latter  case,  it  is  possible  to  define  the  photodesorption  lineshape 
as  a function  of  the  given  light  frequency.  The  desorption  lifetime  can  be  estimated 
from  the  late  time  decay  of  the  autocorrelation  function  — the  projection  of  the 
oesorptron  complex  molecular  wavefunction  on  its  initial  value.  The  lineshape 
is  obtained  as  a temporal  Fourier  transform  of  the  TCF  for  the  transition  dipole 
moment  between  the  electronic  states  of  the  adsorbate-substrate  complex. 

The  desorption  lineshape  expression  in  terms  of  the  dipole-dipole  TCF  is  anal- 
ogous to  the  rate  constant  given  as  the  integrated  collisional  TCF  in  scattering 
theory,  which  was  derived  from  the  Fermi  golden  rule  for  the  transition  proba- 
bility. However,  the  total  energy  in  the  primary  DFs  of  the  adsorbate-substrate 
complex  is  not  conserved  in  this  theory  due  to  the  dissipation;  as  a consequence. 
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■he  transition  probability,  as  well  as  ihe  desorption  lineshape,  acquire  a finite 
width  and  were  demonstrated  to  become  Lorentzian  functions  in  energy  with  the 
peat  shifted  by  the  frequency  of  the  absorbed  light  and  the  width  proportional 
to  the  expectation  value  of  the  dissipative  Hamiltonian.  With  this  generalization, 
the  concept  of  desorption  lineshape  becomes  applicable  to  the  systems  including 
the  energy  dissipation. 

Computational  Algorithm  for  the  CO  Photodesorption  from  the  Ni(OOI)  Surface 

The  CO  photodesorption  dynamics  was  studied  by  solving  the  time  dependent 
Schrddinger  equation  (TDSE)  with  ihe  effective  Hamiltonian  of  the  primary  DFs 
of  the  system.  The  desorption  was  modelled  as  a two  electronic  states'  process, 
with  the  system  evolving  in  two  primary  DFs.  The  two  electronic  stales  were 
described  by  two  complex  potential  energy  surfaces  (PESs)  so  the  effects  of  the 
energy  dissipation  in  die  substrate  were  included.  The  PESs  and  their  couplings 
were  parametrized  semiempirically,  using  the  available  experimental  data  for  the 
Ni(OOI)(ads)CO  system  and  electronic  structure  calculations  for  CO  bound  to  a 
single  N'i  atom  and  an  Ni|0  cluster. 

The  primary  TDSE  was  solved  by  propagation  of  wavepackets.  using  the  split 
operator  propagation  scheme  to  construct  the  time  evolution  operator,  with  a mod- 
ification to  include  the  imaginary  dissipative  potentials.  The  diabatic  electronic 
representation  was  used  in  the  time  propagation,  because  it  eliminates  momentum 
couplings  between  the  electronic  basis  states;  the  real  and  imaginaty  potential 
energy  matrices  were  diagonalized  separately  by  orthogonal  transformations.  The 
time  propagation  algorithm  was  implemented  numerically  using  a discrete  grid 
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lo  represent  the  coordinate  and  momentum  space  — the  variables  of  each  factor 
in  the  split  operator  propagator  were  chosen  such  that  the  Hamiltonian  in  each 
factor  was  diagonal,  and  the  coordinate  and  momentum  space  were  related  by  a 
Fourier  transformation.  A fast  Fourier  transform  (FFT)  algorithm  was  employed 
to  construct  the  kinetic  energy  propagator. 

Summary  of  Results  of  the  CO  Photodcsorption  Calculations 

The  calculations  of  the  CO  photodesorption  in  this  study  were  done  varying 
three  different  parameters:  the  frequency  of  the  absorbed  light  (starting  with  a 
resonant  absoiption),  the  magnitude  of  the  dissipative  imaginary  potentials  and 
their  coupling,  and  the  asymptotic  energy  difference  between  the  two  primary 


because  the  exact  resonance  absorption  frequency,  in  this  case  corresponding  to 
a wavelength  of  240  nm.  may  not  be  experimentally  available.  The  magnitudes 
of  drssipation  potentials  and  their  coupling  were  varied,  as  they  had  been  left 
as  free  parameters  in  the  construction  of  the  model:  neither  previous  references 
nor  the  electronic  structure  calculations  could  provide  their  values.  Finally,  the 
asymptotic  energy  difference  was  varied  because  the  electronic  basis  slates  in  the 
model  represent  the  CO  molecular  orbitals  coupled  with  the  Ni  band  orbitals,  so 
a continuous  variation  of  the  energy  difference  between  those  states  is  possible. 

The  photodcsorption  lineshape  for  the  case  of  resonant  absorption  (fig.  7.6, 
solid  line)  shows  two  peaks  between  4.25  eV  and  5.75  eV;  the  largest  peak  is 
centered  at  the  energy  of  4.50  eV  with  a half  width  of  0.50  eV.  This  is  comparable 
with  the  EELS  spectroscopy  result  [Akimoto  et  al.  1979)  that  shows  a peak  at  5.8 
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eV  wiih  a half  widlh  of  approximately  1 eV.  This  means  that  the  semiempirical 
paramctrization  for  the  PESs  and  their  couplings  used  in  the  calculations  gives 
a good  qualitative  description  of  the  Ni(001)(ads)CO  system,  even  though  the 
potential  energy  functions  were  obtained  with  a single  Ni  atom  or  a relatively 
small  10-atom  duster  simulating  the  solid  nickel. 

The  results  for  the  logarithm  of  the  square  modulus  of  ihe  autocorrelation 
function,  varying  the  absorbed  photon  frequency  (fig.  7.4)  and  the  dissipation 
strength  parameters  (fig.  7.9).  show  that  the  decay  of  the  desorbing  state  is  not 
exponential  at  short  times.  This  is  evidence  that  the  electronic  transitions  in  the 
adsorbate-substrate  complex  occur  both  through  the  diabatic  coupling  of  the  PESs 
and  through  the  dissipative  coupling  that  represents  the  primary  DF  deexcitation 
via  the  energy  dissipation  in  Ihe  substrate.  These  transitions  are  evident  in  the 
results  for  the  relative  population  of  the  ground  electronic  slate  of  the  complex 
(fig.  7.12).  it  is  clear  that  their  probability  is  proportional  to  the  dissipative 
coupling  strength;  even  in  the  absence  of  dissipation  there  is  a finite  probability 
to  find  the  system  in  its  ground  electronic  state  due  to  the  diabatic  coupling  at 
short  adsorbate-substrate  distance. 

From  the  expectation  value  over  time  for  the  adsorbate-substrate  distance  R 
(figs.  7.2  and  7.7).  it  is  possible  to  calculate  the  kinetic  energy  of  the  desorbing 
CO  molecule  and  examine  its  dependence  on  the  absorbed  light  frequency  and 
the  dissipation  strength.  If  the  frequency  of  the  absorbed  photon  is  above  the 
resonance,  the  additional  photon  energy  becomes  the  kinetic  energy  of  the  initial 
state  of  propagation,  with  a momentum 


distribution  that  averages 
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(because  of  the  assumed  Franck-Condon  electronic  transition  upon  absorption). 
As  a consequence,  the  corresponding  wavepacket  is  expected  to  deform  at  short 
times  and  later  evolve  in  the  asymptotic  region  of  large  R with  a larger  velocity. 
The  results  (fig.  7.2)  showed  this  behaviour,  as  the  wavepacket  initially  spends  a 
longer  time  in  the  region  close  to  the  adsorption  site  to  later  move  away  from  it 
faster.  However,  not  all  of  the  extra  energy  from  the  light  absorption  is  converted 
to  the  asymptotic  kinetic  energy  (table  7.II),  due  to  the  presence  of  dissipative 
potentials.  When  the  dissipation  strength  is  increased,  both  the  relative  kinetic 
energy  and  the  internal  vibrational  energy  of  the  desorbing  species  are  reduced 
(figs.  7.7  and  7.8,  table  7.111). 

From  the  results  of  the  square  modulus  of  the  autocorrelation  function  (figs. 
7.4  and  7.9),  it  is  apparent  that  the  desorption  lifetime  evaluated  from  the 
exponcm  I decay  of  the  autocorrelation  function  at  late  limes  bears  no  significant 
dependence  on  the  absorbed  photon  frequency,  at  least  within  the  interval  of  the 
calculations,  however,  it  is  clearly  inverse  proportional  to  the  dissipation  strength 
(table  7.IV).  which  is  the  same  dependence  that  had  been  observed  in  the  model 
calculations  m Chapter  7.  These  results  mean  that  the  energy  of  absorbed  light 
is  distributed  quickly  (within  the  time  at  which  the  adsorbed  molecule  is  still 
close  to  the  substrate)  and.  except  for  a small  change  in  the  kinetic  energy  of 
the  desorbing  species,  it  does  not  significantly  affect  the  late  time  dynamics.  On 
the  other  hand,  the  energy  dissipation  in  the  substrate  plays  a crucial  role  in  the 
photodesorption  dynamics,  both  for  the  lifetime  of  desorbing  states  and  for  the 
internal  vibrations  of  the  desorbing  species  (fig.  7.8). 
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Outlook  for  Future  Research 


This  study  has  developed  a theory  of  photodesotption  that  includes  the 
continuous  energy  spectrum  of  the  substrate  through  the  SCF  coupling,  and 
directly  incorporates  the  energy  dissipation  in  the  substrate  into  the  quantum 
dynamics  of  the  desorbing  species.  It  has  been  shown  that  the  magnitude  of 
dissipation  is  related  to  the  TCFs  of  the  lattice  displacements  and  the  electronic 
density  fluctuations  in  the  substrate  metal,  and  the  computational  results  show  that 
it  significantly  influences  the  motion  of  the  adsorbate.  The  neat  step  should  be 
a detailed  study  of  the  nature  of  the  adsorbate-substrate  coupling  which  should 
give  exact  results  for  the  dissipative  potentials  and  couplings  to  substitute  their 
semiempirical  parametrination  in  this  work.  Also,  the  theory  would  benefit  from 
higher  quality  ab  inilio  calculations  of  the  excited  electronic  stale  PES  lo  improve 
iis  parametrization  and  reproduce  more  closely  the  cxpcrimenlal  results. 

A numerical  algorithm  for  the  wavcpackel  propagation  in  two  degrees  of 
freedom  over  two  coupled  complex  PESs  was  developed  and  applied  on  the  case 
of  CO  desorbing  from  ihe  Ni(001)  surface.  It  allows  for  a quantum  mechanical 
calculation  of  ihe  adsorbate's  motion  relative  to  the  substrate  and  its  internal 
vibration.  However.  Ihe  rotational  motion  of  the  desorbing  species  is  not  taken 
into  account,  so  a further  improvement  of  the  algorithm  would  have  to  include 
another  degree  of  freedom  — a rotation  angle  of  the  diatomic.  A practical  obstacle 
in  this  development  is  the  fact  that  the  computational  time  and  effort  increases 
approximately  quadratically  with  Ihe  number  of  discrete  grid  points  per  each 
degree  of  freedom,  so  it  may  be  necessary  to  employ  more  economical  numerical 


methods,  such  as  the  discrete  variable  representation  (DVR)  [Light  et  al.  1985). 
Also,  the  presence  of  a continuum  of  substrate  electronic  states  may  imply  that 
there  may  be  more  than  one  excited  electronic  state  of  the  adsorbate-substrate 
complex  involved  in  the  dynamics,  so  the  algorithm  should  be  modified  to  include 
another  PES  coupled  diabatically  to  the  existing  ones.  This  would  obviously 
increase  the  size  of  matrices  involved  in  the  computation  and  also  make  the 
transformation  between  the  electronic  representations  more  complicated,  although 
conceptually  unchanged. 

Finally,  it  can  be  expected  that  recent  developments  in  sub-picosecond  laser 
spectroscopy  should  produce  more  high  quality  results  of  time  resolved  experi- 
ments on  the  molecules  adsorbed  on  metal  surfaces.  This  should  provide  more 
information  to  further  improve  the  parametrization  of  model  potentials  and  cou- 


APPENDIX  A 

RATE  CONSTANT  IN  TERMS  OF  A COLLISIONAL  TIME 
CORRELATION  FUNCTION  FOR  A DISSIPATIVE  SYSTEM 


This  appendix  contains  a detailed  derivation  for  the  reaction  rate  constant 
expressed  as  a time  integral  of  a collisional  time  cotTelation  function  (TCF),  for  a 
dissipative  system  interacting  with  an  external  radiation  field.  It  is  a generalization 
of  the  previous  work  (Micha  1981.  Micha  and  Vilallonga  1985]  to  take  into 
account  the  energy  dissipation  in  the  present  system;  in  the  cose  of  interest, 
it  comes  from  the  coupling  of  the  CO  adsorbate  to  an  energy  continuum  of  the 
infinite  Ni  substrate,  and  as  a consequence  the  effective  Hamiltonian  for  the  region 
around  the  adsorbate  contains  an  imaginary  term  and  is  not  a Hermitian  operator. 

One  starts  by  separating  the  system-radiation  term  from  the  rest  of  the  primary 
effective  Hamiltonian 


Ff  = F0  + HpL 
Fo  = ftp  + Hi  — iWp/2 


(A.l) 


where  all  the  operators  have  been  defined  in  Chapter  2.  This  is  similar  to  the 
standard  separation  of  perturbation  theory,  but  in  this  case  there  are  in  principle 
no  limits  to  the  coupling  strength.  The  "uncoupled"  Hamiltonian  Pa  is  obviously 


not  Hermitian,  so  it  defines  a "right"  eigenvalue  problem 


m°„)  = {En-iW„/2M)  (A.2) 


and  a "left"  eigenvalue  problem 


(A.3) 
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where  the  complex  eigenvalues  are  determined  by  £„  = (0<[| (Hr  + 
*t)l#S)/W2l«S>  and  W„  = WS|H,,|^)/(«“|^S).  The  left  and  right  eigenvec- 
tors combine  for  the  orthonormality  and  closure  relations  as  follows 

(#8l*S,)  = *n.m 

EMtioti-/  M 


From  eqns.  (A.2)  and  (A. 3),  it  follows  that  the  time  evolution  oper 
acts  on  the  right  and  left  eigenvectors  producing 

exp  (-iFbl/&)|«“)  = |^S)exp[-i(£„  - iW„/2)t/h] 
(^|exp(iF0«/A)  = exp (>(£„  -iHr„/2)(/ft){^»| 


vith  f0 


(A.5) 


In  analogy  with  Chapter  2,  the  partitioning  (A.l)  of  the  effective  Hamiltonian 
defines  an  interaction  picture  in  which 

WO)/  = exp  (iFot/fi)  |,P(0) 
ft(‘)  = «P  t exp  (-ifbf/ft) 

to  the  primary  TDSE  becomes 


(A.7) 


This  equation  can  be  integrated  formally  to  obtain 


im/  = 1^(0))/  - [ f WWMhdt'  (A .8) 

From  the  definition  of  the  interaction  picture  1^(0)),  = |^(0J).  and  in  the  lowest 
order  of  iteration,  the  interaction  picture  solution  becomes 


!«'))/  = 


jJ  1 1^(0)) 


(A,9) 


For  a particular  choice  of  initial  conditions  |V’P(0))  = |0°).  this  becomes 

NSW)/  = j/  - |#S)  (A.10) 

In  analogy  with  the  case  of  a Hermitian  Hamiltonian,  one  defines  a state  to 
state  rate  constant  [Rodberg  and  Thaler  1967]  as 

Rm~"  = J|«WwW>f  = jflAn—Ml’  (A,  II) 

where  Am—„(i)  denotes  a projection  of  the  interaction  picture  wavcfunction  (A.  1 0) 
onto  the  left  eigenstate  |dS)  of  F0.  Inserting  (A.10)  into  the  equation  for  A„-,„(r), 
using  the  time  dependence  (A.6)  of  F,  and  relations  (A.5)  to  operate  on  the 
eigenstates,  and  finally  integrating  over  time,  the  projection  becomes 


-»(0  = 


(A.  12) 


(E,„  - En)2  + [Wm  - W„f/ 4 
where  the  delta  function  comes  from  the  orthonormaiity  (A.4).  Inserting  (A.  12) 
into  (A.  1 1)  yields  the  state  to  state  rate  constant 


( Wm  - H'„) 


7|  (Em  - &,)'  + - WrfU  (A  I3) 

It  is  noteworthy  that  in  the  absence  of  dissipation,  i.c.  in  the  limit  of  IP,,  and  W„ 
both  going  to  aero,  the  Lorentzian  in  (A.13)  becomes  a Dirac  delta  function  in 
energy,  thus  recovering  the  known  Fermi  golden  rule  result  [Rodberg  and  Thaler 
1967,  Schiff  1968]. 

To  convert  the  rate  constant  (A.  13)  into  an  integral  of  a TCF.  the  Lorentzian 
in  energy  must  be  expressed  as 

- W‘ 


(A.  14) 


‘.JlAIr'M  [nS)«l]f,WI*lW 


APPENDIX  B 

INPUT  DATA  FOR  THE  ZINDO 
CALCULATIONS  OF  NiCO  and  Ni10CO 


This  appendix  contains  information  on  the  input  data  for  the  ZINDO  calcu- 
lations that  were  done  on  the  NitoCO  and  NiCO  systems. 

Specifying  the  atomic  positions  for  CO  adsorbed  on  a Nito  cluster,  it  was 
taken  into  account  that  metal  nickel  forms  a BCC  crystal  lattice;  the  cluster  was 
constructed  from  the  first  and  second  nearest  neighbours  of  the  adsorption  site 
Ni  atom,  and  an  additional  atom  directly  below  the  adsorption  site.  The  atomic 
Cartesian  coordinates  (in  Angstroms)  ate  given  in  table  B.I. 


The  remainder  of  this  appendix  discusses  other  relevant  parameters  for  the 
ZINDO  calculations.  The  SCF  basis  set  for  the  ground  electronic  stale  energy 


calculation  is  specified  according  10  Ihe  atoms  in  the  system;  in  this  case  it  consists 
of  two  atoms  with  s and  p orbitals  (0  and  C)  and  one  (or  10)  atoms  with  s,  p 
and  d orbitals,  for  each  of  the  Ni  atoms.  The  important  distinction  is  made  in 
the  treatment  of  two-electron  integrals:  in  the  Hartrce-Fock  SCF  calculations  for 
the  ground  state  energy  one  must  use  their  theoretical  parameters  ("theoretical 
gammas"),  whereas  in  the  Cl  calculations  of  the  excited  states  the  parameters  that 
are  used  come  from  spectroscopy  ("spectroscopical  gammas").  Other  important 
parameters  are  listed  in  table  B.I1  below. 


Table  B.II:  Input  parameters  for  the  ZINDO  calculations  of  NiCO  and  NiioCO. 


Total  no.  of  valence  electrons 


No.  of  open  shells  in  ROHF-SCF 


APPENDIX  C 

FLOW  CHART  OF  THE  COMPUTER  PROGRAM 
FOR  TIME  EVOLUTION  OVER  TWO  COMPLEX 
PESs  IN  TWO  DEGREES  OF  FREEDOM 
f START  A 
\ program  "cphe2d.f" 
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